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Absorption  Data  Analysis  for  EUVS  Experiment 
on  the  Satellites  AE-C,  D.  and  E 


1.  INTRODUCTION 

The  purpose  of  the  Extreme  Ultraviolet  Spectrophotometer  (EUVS)  absorption 
analysis  system  is  to  determine  quantitative  characteristics  of  the  thermosphere 
(95-400  km)  including  the  diurnal,  seasonal,  and  latitudinal  variations  of  these 
characteristics.  Results  on  particle-number  densities  of  nitrogen  and  oxygen  (N.J( 

O  and  temperature,  total  mass  density,  and  scale  height  have  been  obtained 
for  various  atmospheric  regions.  Other  studies  in  the  future  will  include  the  dis¬ 
tribution  of  vibrational  states  of  molecular  oxygen  and  the  spatial  and  temporal 
variations  of  this  distribution. 

The  EUVS  instruments  take  data  onboard  the  Atmosphere  Explorer  (AE)  satellites 
and  each  consists  of  24  monochromators  which  are  held  in  a  biaxial  girnbal  device 
that  can  be  pointed  at  the  target  position  on  the  solar  disk  commanded  bv  the  experi¬ 
menter.  Twelve  of  these  monochromators,  MN  Nos.  1  through  12,  are  scan-capable, 
that  is,  they  can  scan  a  range  of  wavelengths.  Most  of  these,  that  is,  all  except 
MN  Nos.  11  and  12,  have  full-disk  view  of  the  sun  and  are  primarily  used  for  the 
EUVS  solar  flux  analysis.  The  remaining  12  monochromators,  MN  Nos.  13  through 
24,  are  fixed  at  specific  wavelengths  of  aeronomical  interest  and  have  a  narrower 
field-of-view.  They  are  used  primarily  for  the  absorption  analysis.  A  detailed 
description  was  published  elsewhere.  ' 

(Received  for  publication  8  August  1979) 

1.  Hinteregger,  H.  E. ,  Bedo,  D.  E. ,  and  Manson,  J.  E.  (1973)  The  EUV  spectro¬ 
photometer  on  Atmosphere  Explorer,  Radio  Science  8(No.  41:349-359. 
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2.  K  K.Sl  l.'I'S  OK  M  V  S  \IIS«)KHTIO\  l)\T\  \V\USIS 


Absorption  results  from  El AS  full  into  two  categories.  The  first  covers  re¬ 
sults  based  onl.v  on  observations  from  EUVS.  rhese  are  the  optical  depths.  T  , 
derived  at  each  wavelength  from  the  observed  signals  In  the  relation 


where  S  is  the  observed  flux  and  S  is  the  incident  unattenuated  flux.  These  optical 
depths  arc  obviously  independent  of  any  model  assumptions  on  absorption  cross 
sections  and  provide  the  best  means  of  comparing  any  model  prediction  with  the 
actual  observations  of  ElIVS. 

The  second  category  of  results  are  those  which  arc  based  on  these  optical  depth 
data,  but  do  involve  the  knowledge  of  absorption  cross  sections  and  certain  assump 
turns  of  atmospheric  models,  rhese  results  include  the  aeronomic.il  parameters 
such  as  particle-number  densities  and  temperatures.  While  the  results  of  this 
second  category  may  be  of  more  direct  interest  to  uerononiy,  they  obviously  involve 
the  adoption  of  certain  model  concepts  and  specific  reference  data  required  to 
produce  them.  The  algorithms  which  produce  both  types  of  data  will  be  discussed 
here,  along  with  an  analysis  of  the  types  of  error  encountered  and  the  data  currently 
available. 

Telemetry  data  of  El  VS  observations  come  in  the  form  of  8-sec  blocks  lor  so- 
called  major  frames  of  126  minor  frames  of  1/16  sec  each).  Each  of  these  blocks 
contains  16datagroups  spaced  0.5  sec  apart  ,8  minor  frames'.  The  latter  0.5-sec 
groups  are  called  El  \  S  frames.  The  monochromators  are  assembled  in  pairs 
and  each  pair  (’’module”)  shares  a  common  detector  (either  a  chenneltron  or 
photomultiplier  tube).  This  means  that  only  one  of  the  monochromators  of  a  given 
module  can  be  sampling  solar  radiation  in  any  given  0.  5-sec  interval.  Therefore, 
each  of  the  16  El  \S  frames  of  one  major  frame  contains  not  only  the  count  samples 
from  each  of  the  12  detectors,  but  also  the  aperture-exposure  index  along  with  the 
grating  step  positions  and  other  instrumental  parameters  needed  for  a  complete 
identification  of  the  Eli  VS  data  received. 

Absorption  analysis  deals  with  the  observations  taken  at  JO  selected  wav  elengths 
of  aeronomical  interest.  A  listing  of  those  20  wavelengths  and  the  reference  data  on 
absorption  cross  section  selected  for  the  production  of  density  data  for  the  vehicle 
AE  -l'  are  given  in  Table  1  as  an  example,  where  I  I  1-20  serves  as  the  identification 
number.  The  column  AM!  gives  the  wavelength  (A'  for  each  1  I  number  in  A  units. 
All  20  wavelengths  are  the  same  for  all  three  El  \  instruments  (on  AK-T,  1).  and  1  > 
with  the  exception  of  1.1.  12  which  represents  \  610  A  for  AIM',  but  204  A  for 

A  E  - 1  >  and  E.  Details  of  the  reference  data  will  be  discussed  in  Section  2.  2.  Due 


1 


to  the  aperture  and  step  positions  corresponding  to  the  instrumental  realization  of 
observations  at  these  selected  wavelengths,  only  eight  may  be  sampled  in  any  given 
0.  5-set  interval.  The  overwhelming  part  of  the  existing  EUVS  absorption  data 
reflects  observations  during  experiment  turn-ons  (passes)  with  apertures  left  in  a 
fixed  position.  Therefore  only  8  of  the  20  wavelengths  are  generally  present  in  a 
given  pass.  For  part  of  the  early  life  of  the  AE-C  EUVS  instrument,  the  apertures 
were  alternated  periodically,  and  consequently,  each  1-sec  period  of  two  contiguous 
El  VS  frames  contains  16  wavelengths  in  the  data  for  a  given  pass.  Each  wave¬ 
length  is  then,  of  course,  covered  with  a  time  resolution  of  1.0  sec  (instead  of 
0.  5  sec). 


Table  1.  Twenty  Wavelengths  and  the  Reference  Data  for  AE-C 


1.1 

\NG 

HHLl 

H  itES 

11  UM 

11  LI  M 

mjltlKKt 

lU.2l.OW  Hill 

St  A  1.1  H  I  HI  »• 

M.NMTIX  1 15 

1 

in 

1.5 

460 

2400 

0 

481 

306 

38 

1 

2 

465 

1.5 

50 

800 

0 

527 

352 

47 

1 

3 

584 

0.5 

2100 

6000 

0 

550 

375 

47 

1 

4 

765 

80.  0 

1000 

JO.  000 

0 

500 

330 

38 

l 

5 

855 

50.  0 

400 

9000 

50 

490 

3  20 

47 

1 

6 

1009 

1  28.  0 

3800 

5000 

50 

260 

184 

17 

2 

7 

1457 

1 10.  0 

950 

1400 

50 

280 

204 

18 

1 

8 

1700 

l  20.  0 

6600 

9000 

50 

230 

155 

16 

2 

9 

17  10 

1 10.  0 

9200 

1 1 , 000 

50 

215 

140 

16 

2 

10 

1835 

8.  0 

10. 100 

12.  000 

50 

160 

115 

1  1 

2 

1  1 

256 

2.  0 

80 

500 

0 

520 

345 

4  2 

1 

12 

304 

2.5 

75 

500 

0 

540 

355 

46 

1 

13 

610 

2.  8 

81 

500 

0 

550 

375 

47 

1 

14 

465 

26.  0 

450 

4000 

0 

527 

352 

47 

l 

15 

584 

37.  0 

210 

800 

0 

550 

375 

47 

1 

16 

1026 

25.0 

700 

2600 

10 

235 

159 

15 

2 

t; 

977 

25.0 

600 

2600 

10 

270 

210 

18 

1 

18 

1216 

0.0 

100 

500 

0 

160 

115 

1 1 

2 

19 

1600 

85.  0 

3320 

u000 

50 

260 

184 

18 

2 

20 

1775 

80.  0 

2330 

5000 

50 

180 

125 

12 

2 

Additional  words  in  each  EUVS  frame  are  devoted  to  the  solar  pointing  status 
of  the  EUVS  instrument  and  to  housekeeping  information  such  as  the  high-voltage 
levels.  The  EUVS  absorption  analysis  system  checks  the  solar  pointing  status  words 
and  deletes  from  consideration  major  frames  of  data  taken  when  the  satellite  was 
in  the  dark  or  when  the  sun  was  not  fully  acquired.  It  also  checks  the  point  on  the 
solar  disk  which  is  under  observation  to  ensure  that  the  instrument  is  pointing  where 
commanded.  Pointing  was  directed  to  the  center  of  the  solar  disk  for  all  AE-D  and 
AE-E  geophysical  observations  by  the  EUVS  instruments.  For  AE-C,  central-pointing 
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1/ 


w .1 3  used  only  until  orbit  860.  After  orbit  HtiO  the  pointing  target  was  kept  7  arc 
nun  off  center.  l'his  was  done  to  bring  the  average  alignment  of  the  actual  target 
points  of  the  most  important  monochromators  closer  to  the  true  center  of  the  solar 
disk. 

Another  point  deserves  mention  here.  The  Kl  VS  instrument  on  AK-C  experi¬ 
enced  severe  hardware  problems  with  aperture  masking  and  grating  stepping  which 

2 

have  been  reported  by  flmteregger  et  al.  "  While  these  problems  have  had  a  serious 
effect  upon  the  flux  analysis,  the  absorption  analysis  was  virtually  unaffected,  since 
it  does  not  require  absolute  fluxes. 


t  PKOlil  n  iO!N  OF  OPTICAL  DOTH  D  ATA  (I.MLHO) 

t  l  Inliodlli  Ihiii 

I'he  Kl  \S  instrument  provides  direct  counts  of  the  solar  flux  simultaneously 
for  a  group  of  8  out  of  the  20  fixed  wavelengths  selected  for  analysis.  In  order  to 
prepare  the  raw  data  for  use  in  the  physical  equations,  an  "unattenuated"  reference 
signal  must  be  computed  and  the  count  samples  corrected  for  background  and  filtered 
to  eliminate  aberrant  data  such  as  those  due  to  bad  transmission.  The  data  are 
smoothed  to  improve  the  statistical  significance  and  are  associated  with  a  quantita¬ 
tive  measure  of  data  quality.  I  poll  performance  of  these  functions  a  data  base  of 
the  resulting  optical  depth  data  is  created  and  upgraded  regularly  as  new  telemetry 
becomes  available.  The  resulting  data  files  contain  the  dates  and  times  for  16 
reference  values  of  optical  depth  shown  in  Table  2.  These  dates  and  times  can  be 
used  to  obtain  all  associated  orbit  and  solar  ephemeris  data  needed  to  characterize 
the  absorption  profile  of  the  entire  pass. 


Table  2.  Sixteen  Standard  Optical  Depth  Points  (T  , 


Low 

Attenuation  Region 

Medium 

Attenuation  Region 

High 

Attenuation  Region 

0.  033 

0.  350 

1.  750 

0.  067 

0.  667 

2.  000 

0.  100 

1.  000 

2.  500 

0.  150 

1. 250 

3.  000 

0.  200 

1. 500 

4.  500 

0.  250 

2.  Hinteregger,  H.  K. ,  and  Chaikin,  l..  M.  ( 1977)  EUV  absorption  analysis  of 

thermospheric  structure  from  AP-Satellito  observations  of  1974-1976, 

COSPA  R  Space  Research  17:525-532. 

- —i -  - ww 
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The  first  step  of  data  reduction  is  to  obtain  an  "unattenuated"  reference  signal, 

called  S  . 
o 

The  data  are  scanned  for  an  interval  where  attenuation  is  negligible.  The 
reference  signal  is  then  computed  within  this  interval  and  is  corrected  for  back¬ 
ground.  The  background  corrected  reference  signal  is  then  used  for  all  subsequent 
computations  of  optical  depths.  Once  Sq  has  been  obtained  the  data  are  sequentially 
processed  either  in  the  forward  or  reverse  time  direction.  Each  major  frame  is 
read,  the  signals  corrected  for  background  and  converted  to  optical  depths  by  Eq.  (1), 
that  is, 

T  =  -  tn(S/SQ)  (2) 

where  S  is  the  background- corrected  signal.  These  optical  depth  data  are  fitted  to 
curves  of  the  form  given  by  Eq.  (7).  given  later,  where  £  t  is  the  time  interval  measured 
from  the  start  of  each  fit.  The  range  in  time  over  which  the  curve  fit  is  applied  is 
varied  to  obtain  a  good  signal-to-noise  ratio  and  generally  covers  a  time  period  in 
which  the  remote  probing  point  covers  a  vertical  distance  corresponding  to  a  scale 
height  or  less.  This  regression  provides  the  data  smoothing  required  to  filter  out 
purely  statistical  variability.  The  range  of  the  fit  is  searched  for  the  occurrence 
of  a  standard  optical  depth  value.  If  one  is  found,  the  function  is  inverted,  and  the 
time  corresponding  to  it  is  obtained  and  stored  in  the  data  file.  The  file  also  in¬ 
cludes  the  fitting  coefficients  that  allow  calculation  of  optical  depths  and  the  standard 
error  as  a  measure  of  data  quality.  The  resulting  data  files  can  be  accessed  by  a 
computer  program  (LMC82)  which  associates  the  optical  depths  with  the  correspond¬ 
ing  orbit  and  solar  ephemeris  data. 

3.2  Reference  Signal  Computation  and  background  Correction 

3.  2.  1  THE  SELECTION  OF  "So-EUGIBL.E  KEGIONS" 

The  unattenuated  reference  signal  (Sq)  can  be  obtained  from  the  raw  data  of  the 
satellite  pass  in  several  ways  differing  mainly  in  how  "unattenuated"  is  defined. 

The  values  of  Sq  are  computed  separately  from  each  wavelength,  and  each  wavelength 
may  require  a  different  method  of  computation. 

The  first  and  most  common  method  selects  timespans  for  "regions"  of  the  pass 
so  as  to  meet  certain  criteria  involving  altitude,  solar  zenith  angle,  and  wavelength 
peculiar  absorption  characteristics  included  in  the  reference  data  (see  Table  1). 

For  these  SQ-eligible  regions,  the  reference  signals,  Sq,  are  then  determined  from 
data  within  these  timespans;  the  exact  method  of  this  reduction  is  described  in  the 
next  section.  The  so-called  minimum-ray  height  or  height  of  the  remote  probing 
point  Ov.),  illustrated  in  Figure  1,  is  accessed  once  every  8  sec  (at  the  major  frame 
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start  time)  and  tested  tg.unst  minimum  values  of  altitude  to  determine  if  it  is  within 
the  unattenuated  region.  I’his  height  must  also  be  below  .1  specified  maximum  alti¬ 
tude  to  avoid  signal  contamination  bv  energetic  particles.  This  maximum  is  speci¬ 
fied  in  common  for  all  wavelengths  because  particle  contamination  is  solely  a  func¬ 
tion  of  the  satellite  height.  The  value  of  this  maximum  satellite  height  is  an  input 
parameter  and  can  be  changed  for  each  pass.  In  practice,  the  particle  problem  was 
treated  by  appropriate  selection  of  S  -eligible  regions  at  sufficiently  low  satellite 

heights.  The  minimum  value  of  h. .  allowable  for  determination  of  S  is  different 

w  o 

for  each  wavelength  and  depends  upon  both  the  solar  zenith  angle  and  the  specific 
residual  attenuation  defined  as  negligible  (for  example  T  -  0.  005).  Tor  this 
purpose,  three  reference  parameters  are  associated  with  each  wavelength.  These 
parameters  are  denoted  as  HQH1REF,  llQLOWHET  and  SOU. TUT  HE!  in  Table  1 
corresponding  to  the  notation  of  hy  ;)V  Q.  and  H  respectively,  in  the  text  below. 


Figure  1.  Definitions  of 
Various  Quantities 


The  first  of  the  above  is  the  lowest  acceptable  altitude  for  solar  zenith  angles  greater 
than  05'.  the  second  quantity  is  the  lower  limit  of  altitude  for  overhead  sun  where 
the  column  mass  density  is  a  minimum,  and  the  third  quantity  is  nominal  scale 
height  used  to  calculate  the  altitude  limits  for  any  given  zenith  intermediary  angles, 
lor  angles  less  than  arccosiO.  1),  the  so  called  "sec  \-approximation"  of  the  optical 
depth  factor  is  used,  that  is,  calculating 

\l.XS  Vo*  H  (n  3ec  *  (3) 

where  h  is  the  minimum  of  h  for  the  given  solar  zenith  angle  h, ,  .  the  value 
for  \  =  0  overhead,  and  M  is  the  aforementioned  scale  height  parameter.  The  value 
of  hy  for  zenith  angles  between  arccos(0.  1)  and  05°  is  simply  approximated  by  a 
linear  interpolation  between  the  values  for  \  =  95'  and  the  values  for  \=  arccos(0.  1), 
respectively.  The  specific  values  of  HQHIHET  (h^  s).>,  HQLOWREF  (h^  Q)  and 


SCAL.EHTREF  (H)  shown  in  Table  1  lor  each  wavelength  reflect  the  acceptance  of 
a  residual  optical  depth  of  T  ^  $  0.  005.  The  values  of  these  three  parameters  can 
be  modified  in  accordance  with  any  input  value  of  T  ^  other  then  the  most  commonly 
used  value  of  0.  005.  In  this  case  the  program  calculated  a  corresponding  change 
in  the  acceptable  value  of  h  in  the  form  of 


where  T  „  is  the  selected  value  of  accepted  residual  optical  depth.  Tor  most  of  the 
eccentric  orbital  phases  of  AE-C,  D,  andE,  a  value  of  T  ^  =  0.005  was  used.  How¬ 
ever,  for  the  subsequent  near-circular  orbits,  the  value  of  T  ^  was  raised  to  0.  015 

to  increase  the  rate  of  data  return.  The  actual  value  of  T  is  carried  into  the 

n 

resulting  data  file. 

Depending  on  orbital  characteristics  as  well  as  placement  and  duration  of 
experiment  turn-ons  on  the  AE  satellites,  the  number  of  S^-eligible  regions  ranged 
from  none  to  two  per  pass.  When  two  S^-eligible  regions  do  occur  within  a  given 
turn-on,  the  corresponding  data  pass  was  called  "double  pass.  "  The  data  reduction 
system  can  seek  either  one  or  two  S  -eligible  regions  for  each  wavelength  sepa¬ 
rately,  or  the  experimenter  can  manually  enter  the  regions  defined  as  unattenuated 
for  any  selected  wavelengths.  In  the  latter  case,  the  system  will  then  accept  these 
regions  without  executing  the  aforementioned  calculation  of  a  minimum  acceptable 

value  of  h„. 

W 

When  AE-C  and  AE-E  went  into  circular  orbits  at  relatively  low  altitude 
(-  250  km),  data  on  many  of  the  strongly  attenuated  wavelengths  were  lost  due  to  the 
lack  of  truly  unattenuated  SQ-eligible  regions.  To  improve  the  rate  of  data  return 
over  timespans  of  aeronomical  interest,  the  method  of  direct  input  of  the  timespan 
for  determination  of  Sq  then  had  to  be  used  frequently.  The  SQ  values  so  derived 
were  of  course  not  truly  unattenuated  but  reflect  actual  optical  depths  ranging  up  to 
0.  3  for  the  w-orst  cases.  These  values  of  S  could  be  corrected  later  by  a  different 
section  of  the  data  reduction  system  by  calculating  the  residual  optical  depth  for  the 
time  of  the  SQ  observations,  T  (SQ),  to  a  degree  of  accuracy  considered  quite 
adequate  in  view  of  the  relatively  small  magnitudes  of  T  (SQ)  for  most  wavelengths. 

In  some  rare  instances,  the  experimenter  entered  actual  Sq  values  directly 
instead  of  merely  selecting  the  timespan  for  which  S  would  then  be  determined  by 
the  program. 

Regardless  of  how  timespans  of  So-eligible  regions  are  determined,  the  actual 
existing  measurements  are  screened  for  insufficient  data.  For  instance,  regions 
having  an  insufficient  amount  of  TM  data  are  rejected.  This  rejection  threshold  is 
an  input  parameter  that  was  normally  set  to  require  a  minimum  of  one  major  frame 
(8  sec).  It  should  be  noted  that  S  values  are  not  taken  from  the  whole  S  -eligible 


region  whenever  the  latter  is  very  large.  The  maximum  size  of  that  part  of 
So-eligible  regions  for  which  is  actually  calculated,  represents  another  input 
parameter  (routinely  left  at  five  major  frames  for  data  production).  When  it  was 
desired  to  avoid  significant  backgrounds  by  energetic  particle  fluxes  during  the 
period  of  satellite  orbits  with  high  apogee  altitudes,  the  actually  used  sections  of 
large  .iigible  regions  was  selected  to  be  as  close  to  earth  as  possible.  Alter¬ 
natively,  where  particle  contamination  is  not  a  problem,  for  example,  during  the 
phase  of  circular  orbits  with  altitudes  below  400  km,  the  instructions  were  switched 
to  s?lect  the  highest  part  of  the  S^-eligible  region.  In  the  routine  data  production, 
each  of  these  two  options  was  used  where  it  would  give  the  best  results,  an  identifica¬ 
tion  of  this  option  was  carried  into  the  resulting  data  file. 

3.  2.  2  THE  HAW  SIGNALS  IN  S  -ELIGIBLE  REGIONS 

o 

All  of  the  data  points  within  the  selected  section  of  an  S^-eligible  region  are 

input  to  a  standard  least  square  linear  regression.  The  independent  variable  is  the 

time  elapsed  in  seconds  (measured  from  the  start  of  the  selected  section).  Cramer's 

rule  has  been  used  to  solve  the  resulting  two  equations  for  two  unknowns.  Singular 

values  were  identified  by  a  comparison  of  the  absolute  value  of  the  determinant  with 

-  1 5 

a  minimum  value.  This  tolerance  is  an  input  parameter  which  was  set  at  t.  OX  10 
for  routine  data  production.  Additionally,  another  input  parameter  defined  the 
minimum  number  of  points  for  which  an  SQ  determination  is  considered  successful. 
This  parameter  was  set  at  six.  Data  fits  having  too  few  points  or  singular  deter¬ 
minants  were  rejected,  and  further  processing  for  that  wavelength  was  then  omitted. 

For  successful  regressions,  the  fitted  value  of  the  raw  signal,  Ho(Sq),  was  then 
extracted  for  the  center  of  the  fitting  interval.  The  change  in  signal  over  the  whole 
S^-fit  range  is  a  measure  of  the  quality  of  SQ  and  is,  therefore,  included  in  the  data 
file  records.  Ideally,  it  should  be  zero. 

To  protect  the  results  from  contamination  by  abberrant  data,  some  form  ^f 
data  filtering  for  the  raw  signal  computation  is  obviously  required.  This  is  accom¬ 
plished  by  comparison  with  fixed  upper  and  lower  limits  given  as  reference  data  for 
each  wavelength.  These  are  set  to  take  a  wide  range  of  values,  and  are  essentially 
designed  to  filter  out  transmission  errors,  stuck  bit  errors,  and  so  on.  The  values 
for  these  limits  are  in  Table  1  (see  quantities  called  lULIIVl  and  1LLIM).  Once  the 
raw  signal  (RQ)  is  computed,  the  background  is  subtracted  and  the  unattenuated 
reference  signal  is  given  as  SQ  =  Rq  -  B. 

3.  2.  3  THE  BACKGROUND  CORRECTION 

A  correction  for  scattered-light  background  is  accomplished  by  the  subtraction 

S  =  R  -  B,  where  B  is  determined  by  a  reference  value  of  B,,,,t,  stated  for  a 
o  o  J  RLr 

reference  value  of  SQ  =  SQ  ^p.  *’°r  a  given  value  of  Rq,  the  corresponding  back¬ 

ground  is  scaled  appropriately,  that  is. 
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(5) 


QREF  Ro/(So,  ref  *  ‘’kef* 


where  B,,,.,.  and  S  are  listed  for  each  wavelength  in  Table  1  (denoted  as 

REr  o,  REr 

BREF  and  SOREF  respectively).  Changes  in  instrumental  sensitivies  were  caused 
by  the  changing  detector  responses  and  occasional  stepping-up  of  the  operational 
high  voltage.  The  latter  was  undertaken  to  compensate  for  degradation  of  pulse- 
height  distributions.  Since  the  optical  depth  data  are  generated  by  taking  ratios  of 
the  signals,  r  =  fn(S  /8),  the  absorption  analysis  is  obviously  not  sensitive  to 
long-term  variation  of  the  absolute  photometric  sensitivities. 

3.  2.  4  FINAL  COMPUTATION  OF  SQ 

The  final  value  of  the  unattenuated  reference  signal  is  computed  by  subtracting 
the  scaled  background  from  the  computed  raw  signal  as  shown  in  Eq.  (5).  For 
single  passes  (satellite  data  passes  containing  only  one  So-eligible  region)  this  is 
the  S  value  used  for  all  optical  depth  computations  for  that  pass.  In  double  passes, 
where  two  slightly  different  reference  signals  are  obtained  for  two  parts  of  the  pass, 
the  final  S  used  for  the  optical  depth  computation  is  obtained  by  a  linear  interpola¬ 
tion  between  the  two  values  for  the  point  of  time  measuring  the  attenuated  flux  S. 

The  magnitude  of  the  unattenuated  average  count  sample,  R^IS^),  effectively 
determines  the  upper  limit  of  optical  depths  beyond  which  no  significant  data  can  be 
expected.  The  precision  to  which  optical  depths  may  be  determined  is,  of  course, 
also  related  to  the  magnitude  of  Rq.  In  the  interest  of  data  quality,  the  system 
therefore  does  not  produce  optical  depth  data  whenever  the  background-corrected 
signal  becomes  less  than  five  counts  per  0.5-sec  sample.  This  means  that  the 
system  will  seek  optical  depths  only  up  to  a  certain  highest  applicable  standard 
optical  depth  reference  value,  that  is. 


T  ....  <  fn  (S  / 5)  .  (6) 

\1  A  A  o 

Hence  this  maximum  is  different  for  different  wavelengths.  If  S  for  any  wavelength 
is  so  low  that  none  of  the  standard  optical  depths  will  be  sought,  processing  for  that 
wavelength  is  halted.  Similarly,  processing  is  also  terminated  for  wavelengths 
whose  S^-eligible  region  extends  over  the  whole  pass. 

3.3  Opliral  Depth  Lotiipiinilion 

3.  3.  1  GENERAL  DESCRIPTION 

Optical  depth  data  in  the  absorption  analysis  are  given  in  terms  of  the  times  of 
occurrence  of  optical  depths  corresponding  to  16  specific  reference  values.  These 
reference  values  are  given  in  Table  2.  Orbital  parameters  may  be  computed  for 
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each  time,  and  so  the  optical  depth,  altitude,  latitude,  longitude,  and  so  on,  arc  all 
associated  via  this  time.  Scale  heights  and  time  variations  mas  be  inferred  from 
differences  in  these  data,  but  such  computations  will  not  be  discussed  here.  The 
most  important  reason  for  the  use  of  a  standard  optical  depth  list  is  the  substantial 
reduction  of  data  quantity.  Using  this  method,  the  information  from  an  entire  pass 
is  compressed  into  data  on  Iti  points  while  losing  very  little  information.  This 
enables  follow-on  processing  to  proceed  without  requiring  excessive  amounts  of 
computer  time  or  storage.  In  fact,  the  ability  of  the  KUYS  absorption  system  to 
scan  through  large  time  periods  of  existing  measurements  is  due  to  this  data  com- 
pa  rison. 

The  KUYS  data  reduction  system  starts  with  either  the  smallest  or  greatest  of 
the  standard  optical  depth  values  depending  on  whether  the  actual  attenuation  is 
increasing  or  decreasing.  1'he  attenuation  trend  is  known  because  the  data  reduction 
system  finds  the  time  corresponding  to  the  lowest  probing  height  ill  ,1  in  the  pass 

w 

(tli  tt  is,  close  to  point  of  highest  attenuation)  and  has  either  computed  or  has  had 
as  input  the  direction  of  processing.  In  order  to  seek  the  largest  or  smallest 
optical  depth  point,  the  system  performs  a  two-parameter,  non-linear  exponential 
regression  of  the  form 


T  A  ( 


UAt 


(7) 


where  A  and  It  are  the  fitting  coefficients,  and  At  is  the  time  difference  in  seconds 
from  the  starting  point  of  the  regression.  The  time  inverval  over  which  this  re¬ 
gression  is  taken  depends  upon  the  amount  of  signal  change.  Tor  periods  of  rapid 
change,  the  interval  is  contracted,  so  that  important  features  are  not  lost;  for 
periods  of  slower  change,  the  interval  is  expanded  to  improve  the  statistical  accu¬ 
racy.  The  details  of  determining  this  interval  arc  given  in  the  next  section.  The 
T  values  are  computed  from  the  observed  signals  by  Kq.  (2),  where  S  is  the  back¬ 
ground-corrected  signal.  When  the  data  fit  lias  been  completed,  the  endpoints  tin 
optical  depth)  are  checked  to  see  if  the  reference  point  is  contained  in  the  interval. 

If  so,  tlu>  relation  of  Kq.  (2)  is  inverted,  and  At  is  computed  for  the  standard  optical 
depth  point  being  sought.  Addition  of  At  to  the  starting  time  of  the  interval  then 
produces  the  time  corresponding  to  the  occurrence  of  one  of  the  standard  values  of 
optical  depths.  The  next  optical  depth  point  in  the  list  then  becomes  the  one  sought, 
and  tin'  current  interval  is  checked  to  see  if  that  point  is  within  the  present  fitting 
range.  If  not,  the  fitting  interval  is  shifted,  another  regression  performed,  and  so 
on,  until  the  next  point  is  reached.  The  fitting  intervals  overlap  for  part  of  their 
ranges  to  insure  a  smooth  transition  between  regressions  and  the  degree  of  overlap 
depending  on  the  time  interval,  l'lio  interval  lengths  (overlap  times)  are  giv  en  as 
one  of  the  three  sets  24(8)  sec,  16(8),  8(4)  sec  where  the  numbers  in  parentheses 
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expose  the  overlap  of  successive  fit  intervals.  When  optical  depth  points  are  not 
attained  or  are  skipped,  such  as  when  the  satellite  is  always  too  high  (or  low)  for 
the  hy  to  correspond  to  those  points,  the  system  notes  the  condition  and  proceeds. 

As  part  of  the  results  for  each  standard  point,  the  standard  error  is  computed  and 
preserved  as  a  measure  of  data  quality.  The  fitting  coefficients  are  also  kept.  For 
double  passes,  two  records  are  produced,  one  for  the  upleg  and  one  for  the  downleg, 
and  for  each  leg  all  16  standard  points  may  be  computed. 

3.  3.  2  RANGE  SELECTION  FOH  OPTICAL  DEPTH  REGRESSION 

When  the  satellite  is  over  the  sunward  side  of  the  earth  (solar  zenith  angle 
-  90°),  the  minimum  ray  height,  h,.,  is  falling  or  rising  only  as  fast  as  the  satellite 
itself,  and  so  the  changes  in  signal  are  relatively  slow.  However,  as  the  satellite 
goes  into  or  out  of  eclipse,  the  minimum  ray  height  may  fall  or  rise  very  rapidly. 
Obviously  the  rate  of  change  is  most  rapid  when  tfie  sun  is  close  to  the  plane  of  the 
satellite  orbit.  The  situation  is  then  similar  to  the  data  acquisition  in  a  rocket 
experiment.  To  prevent  statistical  fluctuations  from  distorting  the  results  in  the 
first  case,  and  to  prevent  the  loss  of  information  caused  by  the  use  of  the  exponential 
in  the  second  case,  different  sizes  of  fitting  interval  are  required.  The  exponential 
function  was  used,  rather  than  a  linear  or  polynomial,  because  it  more  closely 
represents  the  true  physical  function,  and  thus  can  cover  a  wider  range.  This  func¬ 
tion  represents  the  data  well  over  ranges  where  the  change  in  signal  corresponds  to 
nearly  a  full  height  of  penetration  into  the  atmosphere.  To  meet  this  requirement 
the  fitting  interval  is  shortened  as  the  level  of  attenuation  increases,  because  the 
scale  height  is  then  also  decreasing. 

The  speed  at  which  the  signal  changes  depends  on  three  factors:  the  wavelength; 
the  current  level  of  attenuation;  the  range  of  solar  zenith  angles  ("x- region"). 

The  wavelength  dependence  is  caused  by  the  different  values  of  absorption  cross 
sections  listed  in  Table  3.  For  this  analysis  it  is  sufficient  to  divide  the  wavelengths 
into  two  groups,  that  is,  more  rapidly  attenuating  wavelengths  indicated  by 
SENSITIVITY  =  2  in  Table  1  and  less  rapidly  attenuating  wavelengths  indicated  by 
SENSITIVITY  =  1  in  Table  1.  The  corresponding  distinction  is  part  of  the  reference 
data  set  for  that  wavelength. 

The  type  of  attenuation  region  (low,  medium,  high)  for  the  standard  optical 
depths  being  sought  is  the  second  factor  determining  the  length  of  the  fitting  range. 

At  higher  optical  depths,  the  scale  height  for  a  given  wavelength,  H(T^),  is  generally 
smaller  than  the  value  H(r^)  for  lower  optical  depth  (lower  kinetic  gas  temperature 
for  lower  altitudes).  Therefore,  the  list  of  standard  optical  depths  is  divided  into 
three  sections:  Low  attenuation  region  (6  values  of  r  ^p,  optical  depth 
0.  033  -  T  5  0.  25),  medium  attenuation  region  (5  values  of  t  ^p,  optical  depth 
0.  35  i  T  5  1.  50),  and  high  attenuation  region  (5  values  of  r  optical  depth 

1.75  s  T  5  4.  5)  as  shown  in  Table  2. 
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Table  .1.  Absorption  Cross  Sections  for  14  of  the  20  Standard  Wavelengths  for 
Which  Aeronomical  Kcsults  Have  been  Created  So  Par 


LL 

\ 

nom 

(A) 

C  ross 

O 

Sections  (Mb) 

°2  N2 

Hemarks  or  note  of  main  absorbers 

1 

175 

3.5 

5.  1 

O  and  N2 

2 

465 

10.  2 

21.8 

Much  poorer  than  for  No.  14 

3 

584 

13.  0 

23.  0 

O  and  N., 

4 

765 

Problems  of  radiation  hardening 

5 

855 

7.  2 

13.4 

36.  0 

O  and  N^ 

6 

1609 

4.  7 

°2 

7 

1457 

14.  0 

°2 

8 

1700 

(),,  vibr.  distrib.  study 

9 

1730 

(.)„  vibr.  distrib.  study 

10 

1835 

O,.  band  abs.  (difficult) 

1 1 

256 

7.  1 

9.  6 

O  and  N^ 

12 

304 

9.  0 

1 1.6 

O  and  Nr) 

13 

610 

13.  1 

23.  1 

O  and  N., 

14 

465 

10.  2 

21.  8 

O  and  N0 

15 

585 

13.  0 

23.  0 

Poorer  than  for  No.  3 

16 

1026 

1.5 

0.,  (see  Note  1) 

17 

977 

Problems  of  radiation  hardening 

18 

1216 

0.01 

C >„  (see  Note  2) 

19 

1600 

5.7 

(most  observations) 

20 

1775 

(>2  vibr.  distrib.  study 

* 

Note  1:  Absorption  of  11  l,y-/l  at  higher  altitudes  is  complicated  by  absorption 
duo  to  atomic  oxygen. 

Note  2:  The  use  of  the  nominal,  constant  value  of  0.01  Mb  is  justified  as 
first  approximation  only. 

The  final  factor  taken  into  account  is  the  solar  zenith  angle.  As  was  previously 
mentioned,  the  rate  of  change  of  r  increases  with  the  solar  zenith  angle  \.  Again, 
the  range  of  \  is  split  into  three  regions: 

Xji  Xs  arccos  (0.  1 ) 

X2:  arccos  (0.  1)  <  x  <95° 

X3=  X  *  95  • 
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The  type  of  range  chosen  for  :> ny  rog rosston  depends  upon  the  throe  factors 
mentioned  curlier.  Although  the  experimenter  could  modify  the  column 
Sl-.NSl  m  ITV"  of  l  ublc  I,  in  all  routine  data  I'cductvons  the  fitting  intervals  were 
fixed  as  given  in  Table  4.  As  the  lowest  optical  depths  in  the  low  attenuation  region 
are  reached,  say  T  i  0,  1,  the  34  see  width  of  the  widest  range  can  become  Inade¬ 
quate.  I  or  this  reason,  the  two  lowest  standard  optical  depth  points  are  often 
missing.  However,  if  the  count  sample  for  the  specific  wavelength  is  high  t  reducing 
statistical  scatter)  good  values  may  he  obtained  even  forth*’  lowest  standard  value 
of  r  |j|,  0.03d.  With  the  exception  of  the  lowest  t  values,  the  method  of 

range  selection  does  keep  the  statistical  variability  and  true  functional  shape  within 
an  appropriate  cause  for  the  regressions.  Additionally,  the  number  of  data  points 
included  within  each  fit  is  kept  high  enough  t * >  avoid  spurious,  physically  lion- 
meaningful  solutions  and  yet  low  enough  to  facilitate  rapid  processing  of  the  several 
hundred  regressions  which  are  performed  fur  each  pass. 

To  provide  as  much  uniformity  to  the  regression  as  possible,  three  standard 
fitting  types  have  been  selected.  They  are  listed  in  Table  -l  under  various  com¬ 
binations  of  tin'  three  factors  discussed  above.  I'll*'  first  designated  34i!'l  above, 
uses  a  range  of  three  major  frames  (34  see.  up  to  It'  data  points)  with  one  major 
frame  overlap  it'  sec.  16  points)  of  the  previous  regression  range;  the  second, 
denoted  16(6),  uses  two  major  frames  (16  sec,  3c  points)  one  of  which  overlaps  the 
previous  range;  the  third,  designated  )'i4).  uses  one  major  frame  it'  see.  16  points) 
and  one-halt  of  a  major  frame  as  overlap  14  see,  ti  points).  The  number  of  points 
mentioned  above  refer  to  the  maximum  possible  that  will  occur  only  for  perfect 
I'M  data  and  constant  aperture  selection.  When  the  Kl'VS  instrument  used  aperture 
alternation  (AK-C  orbits  up  to  orbit  6601  16  wavelengths  instead  of  6  wore  sampled 
and  Hi*'  corresponding  numbers  of  points  per  wavelength  are  then  halved.  Most 
often,  the  maximum  is  indeed  attained,  because  the  occurrence  of  bad  telemetry 
over  a  time  interval  as  short  as  34  sec  would  have  to  be  very  severe  to  have  a  major 
effect  on  more  than  one  wavelength.  Fortunately,  this  has  rarely  happened. 

3.3.3  UKU  MISSION  ALGOMTHM 

Tin'  algorithm  used  to  fit  Kq.  (7)  was  a  standard  non-linear  least  squares  fit 
to  th*’  exponential.  A  linear  fit  to  the  logarithm  of  optical  depth  could  not  be  used 
because  the  optical  depth  values  could  be  negative.  While  a  negative  optleal  depth 
has  no  physical  significance  (the  atmosphere  does  not  add  ultraviolet  radiation), 
such  points  can  occur  because  of  normal  statistical  variability  of  signals  in  the  low- 
attenuation  region.  To  reject  these  points  would  (hus  cut  off  statistical  fluctuations 
in  one  direction  while  leaving  them  intact  in  th*'  other,  thus  distorting  the  statistical 
treatment  of  the  data. 


Non-linear  regressions  require  an  initial  approximation.  Like  many  other 
functions,  the  one  used,  required  the  approximation  to  be  reasonably  close  to  the 
result.  The  data  points  are  collected  monotonicully  in  time  (and  thus  in  attenuation), 
and  the  initial  approximation  is  that  exponential  which  would  generate  the  final 
optical  depth  in  the  range  from  ttie  initial  ones.  Thus,  the  approximation  for  A  is 
the  average  of  the  first  three  optical  depth  values  to  represent  the  initial  point,  and 
the  approximation  for  B  is  the  value  required  to  generate  the  last  point,  which  is 
represented  by  the  average  of  the  final  three  values.  Since 

Average  (final)  =  Average  (initial)  exp  [  IV  (RANGE)]  , 

B  =  tn  [  (Average  (final)/Average  (initial))]  /  Range  . 

This  approximation  procedure  is  very  good  in  all  regions  having  attenuated  data, 
and  speeds  the  regression  considerably.  However,  in  areas  that  are  in  the  very 
low  attenuation  region,  it  sometimes  occurs  that  the  approximation  for  A  would  be 
negative.  Since  ttiis  would  not  lead  to  a  good  solution,  the  approximation  of 
A  =  0.01.  B  =  0.04  is  used  whenever  this  occurs. 

The  iterative  process  of  a  non-linear  regression  requires  a  number  of  statisti¬ 
cal  controls  to  ensure  tin-  validity  of  the  results.  Data  filtering  for  erroneous  data 
will  be  treated  in  Section  3.  3.4,  while  the  control  of  the  fitting  algorithm  and  the 
quality  checking  of  its  results  will  be  covered  here.  To  ensure  against  a  singular 
system  of  equations  or  an  insufficient  number  of  points,  the  determinants  computed 
and  the  number  of  data  points  within  each  regression  are  checked  against  preset 
minimum  values  which  are  the  s  me  as  those  used  for  the  reference  signal  (S  ) 

determination.  The  tolerance  for  ending  the  iterative  process  is  an  input  parameter 

_  £ 

which  for  routine  data  production  was  set  at  IX  10  '  .  The  maximum  number  of 
iterations  permitted  was  20.  If  any  of  the  above  conditions  were  not  satisfied,  the 
fit  was  rejected.  Upon  completion  of  the  regression,  the  resulting  coefficients  were 
compared  with  the  initial  approximation  values.  If  either  of  the  coefficients  was  less 
than  the  initial  approximation  bv  a  factor  of  0.  01,  or  exceeded  it  by  more  than  a 
factor  of  100,  the  fitting  results  were  rejected.  This  was  done  to  prevent  the 
occurrence  of  physically  non- meaningful  solutions.  The  factor  of  100  was  chosen, 
because  in  the  low  attenuation  region,  fairly  large  factors  are  reasonable.  In  addi¬ 
tion  to  these  tests,  another  feature  of  quality  assurance  of  the  system  deserves 
mention  here.  If  the  reference  signal  is  low,  the  higher  optical  depths  would  corres¬ 
pond  to  a  very  small  number  of  counts,  leading  to  imprecision.  For  this  reason,  the 
system  does  not  produce  data  for  the  highest  values  of  r  for  certain  wavelengths 
for  which  the  result  would  be  based  on  an  essentially  insignificant  number  of  count 
samples.  The  method  of  achieving  this  was  covered  in  Section  3.  2.4.  Within  the 
frame  of  quality  checks  and  data  filtering  discussed  next,  regions  of  relatively  poor 


data  quality  can  still  be  used  to  extract  whatever  useful  information  may  be  con¬ 
tained,  without  risking  inadvertent  acceptance  of  erroneous  results  as  a  consequence. 

For  values  of  t  for  which  data  are  produced,  the  latter  include  identification 

HEr 

of  the  standard  error  computed,  to  provide  quantitative  information  on  data  quality. 

3.3.4  DATA  FILTERING 

Data  filtering  for  transmission  errors,  stuck  bits,  and  soon,  is  more  complex 
for  the  optical  depth  computation  than  for  the  reference  signal  determination.  As 
the  count  sample  drops,  the  magnitude  of  the  error  that  will  damage  the  analysis 
also  drops.  This  means  that  for  high  signals,  a  low-order  bit  error  will  not  affect 
the  results  significantly,  that  is,  the  effect  will  be  no  more  than  that  of  the  normal 
statistical  fluctuation  of  the  error-free  data.  However,  as  the  signal  drops  due  to 
attenuation,  the  same  low-order  bit  error  may  correspond  to  an  erroneous  signal 
which  could  be  several  orders  of  magnitude  greater  than  the  actual  signal.  There¬ 
fore,  as  the  signal  drops,  the  risk  of  significant  errors  increases  and  it  becomes 
necessary  to  have  a  dynamic  set  of  limit  filters.  This  is  discussed  next. 

The  absorption  analysis  data  reduction  system  uses  a  pair  of  filters  for  each 
wavelength,  one  high  limit  and  one  low  limit.  Each  limit  is  essentially  a  predicted 
signal,  plus  or  minus  several  deviations,  and  is  changed  whenever  the  regression 
has  accomplished  data  for  one  of  the  standard  optical  depth  points.  The  signal  pre¬ 
diction  for  these  limit  computations  is  subject  to  the  same  factors  that  determine 
the  fitting  interval,  that  is,  the  wavelength,  the  attenuation  level,  and  the  rate  of 
change  of  the  latter.  The  distinction  of  wavelengths  has  already  been  taken  into 
account  because  each  wavelength  has  its  own  set  of  upper  and  lower  limit  filters. 

Since  the  optical  depth  is  a  measure  of  the  attenuation  level,  all  that  is  required  for 
the  optical  deptli  prediction  is  an  estimate  of  rate  of  change  of  the  optical  depth. 

This  estimate  is  based  on  the  solar  zenith  angle.  For  not  too  large  zenith  angles, 

X  <  arccos  0.  1,  the  adopted  system  assumes  that  the  signal  may  diminish  to  a  value 
as  low  as  that  for  three  steps  of  r  above  the  value  of  r  just  attained.  That 
is,  when  one  optical  depth  point  has  been  filled,  the  next  is  then  sought  and  the  new 
low  limit  is  adjusted  to  pass  appropriately  lowered  signals.  The  corresponding 
optical  depth  displacement  for  the  ranges  arccos  0.  1  <  \  v  5)5°  and  \  >  95°,  is  four 
and  five  steps,  respectively.  The  predictions  need  never  be  set  to  an  optical  depth 
higher  than  the  maximum  of  T  =4.5  listed  (see  Table  2).  When  that  point  is 
filled,  no  more  data  of  that  wavelength  in  the  given  pass  will  be  processed.  The 
actual  filter  progression  is  based  on  passing  variations  within  a  stated  permitted 
range  of  deviations  as  discussed  below. 

Once  the  estimate  of  projected  value,  r  is  obtained,  it  is  used  to  determine 
a  signal  by  the  equation 

S  ,  =  S  e  T  est  .  (8) 

est  o 
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F rom  this  signal,  the  lower  limit  is  obtained  by  subtracting  a  certain  number  of 
deviations.  Each  deviation  is  defined  as  the  square  root  of  the  predicted  low  signal. 
The  actual  number  of  deviations  is  an  input  parameter  but  was  set  equal  to  5  for  all 
routine  data  production. 

The  estimate  of  an  appropriate  upper  limit  is  obtained  by  using  the  T  ^,^,-point 
just  attained  (with  no  offset).  The  new  upper  limit  is  then  obtained  by  adding  the 
aforementioned  number  of  deviations  to  the  signal  determined  for  the  last  r 
value.  The  process  for  updating  limits  when  the  optical  depth  is  decreasing  (sunrise) 
is  analogous.  However,  now  the  upper  limit  has  the  offsets  of  three,  four,  or  five 
steps  in  r  and  the  lower  limit  is  based  on  the  T  ^gp-point  just  attained.  Again, 
it  is  not  necessary  to  go  beyond  the  boundaries  of  the  T  jjgj.-  list  (see  Table  2).  If 
at  any  time  the  predicted  lower  limit  of  a  signal  is  less  than  60  counts,  that  limit 
is  set  to  zero.  Immediately  after  the  computation  of  the  unattenuated  reference 
signal,  the  limits  are  set  as  wide  as  possible  (lower  is  zero,  upper  is  So  plus 
5  deviations)  because  the  first  optical  depth  point  to  be  attained  is  unknown  until  it 
is  reached  and  the  limit  updating  procedure  above  takes  over.  Each  time  a  violation 
of  one  of  the  limits  occurs,  the  point  is  rejected  and  the  recorded  count  of  violations 
is  increased  by  one.  Separate  counts  are  maintained  for  upper  and  lower  limit 
violations.  There  is  also  a  counter  for  "fill"  data,  that  is,  where  TM  data  of  geo¬ 
physical  significance  do  not  exist  (having  been  replaced  by  a  specific  flag  indicating 
that  condition).  All  three  of  these  counter  values  are  carried  into  the  data  file 
record  produced. 


4.  PRODUCTION  OF  ALTERNATIVE  OPTICAL  DEPTH  DATA  (I.MC85) 


4.1  liitmtlurtion 


The  EUVS  absorption-data  production  program  named  LMC85  creates  records 
of  considerably  greater  detail  than  those  produced  by  the  standard  program  named 
LMC80  which  creates  data  for  up  to  sixteen  standard  values  of  optical  depths  only. 
The  two  programs  differ  in  several  important  aspects.  The  program  LMC'85  creates 
time-ordered  records  of  optical  depths  for  input-controlled  fixed  stops  in  time, 
whereas  LMC'80  creates  records  for  the  times  of  occurrence  of  certain  standard 
optical  depths  only.  Typical  values  of  the  fixed-time  increment  chosen  for  LMC85 
are  5  or  15  seconds,  selected  to  facilitate  the  comparison  with  the  so-called 
"United  Abstract  Data"  (UA  data)  on  results  of  other  AE-satellite  experiments  given 
with  a  standard  granularity  of  15  seconds.  Results  of  EUVS  observations  are  not 
conducive  to  this  type  of  presentation,  mainly  because  of  the  need  for  a  total  number 
of  words  per  tlmepoint  in  excess  of  the  acceptable  allotment  for  each  experiment. 


Therefore,  it  was  decided  that  standard  UA  file  records  would  not  be  written  for 
the  EUVS  experiment. 

Therefore,  the  records  created  by  the  program  LMC85  represent  a  superior 
substitute  for  otherwise  non-existent  UA  file  data  on  EUVS  absorption  results. 
However,  for  economical  reasons,  no  attempt  can  be  made  at  any  quasi-duplicate 
application  of  I.MC85  approaching  even  remotely  the  complete  coverage  provided 
by  the  standard  optical-depth-data  production  of  LMC80.  Fortunately,  the  latter 
represents  quite  adequate  scientific  information  for  most  purposes.  Hence,  the 
execution  of  the  program  1.MC85  is  expected  to  remain  restricted  to  special  cases. 

The  output  of  LMC85  represents  "stand-alone  files.  "  This  means  that  these  file 
records,  written  in  EBCDIC  format,  contain  all  of  the  important  geophysical  param¬ 
eters  of  interest  for  each  point  in  time  for  which  the  optical  depth  was  evaluated, 
that  is.  the  effort  of  scientific  access  to  these  results  (stored  on  magnetic  tape) 
does  not  require  any  read  program  with  access  to  any  additional  data  such  as  orbit- 
attitude  or  solar  ephemeris  data.  The  format  of  these  files  is  specified  in  Section  4.  3. 

Unlike  LMC80,  the  program  LMC85  produces  optical  depth  data  for  uniformly 
spaced  time  points  of  the  entire  period  of  observations  made  during  a  given  pass. 

Both  the  unattenuated  reference  signal  (Sq)  and  the  optical  depths  are  computed 
from  time-based  linear  regressions.  The  optical  depth  data  reduction  can  be 
specified  to  start  at  some  specific  data  and  time,  with  optical  depths  being  extracted 
for  uniform  time  intervals  which  are  usually  specified  as  15  or  5  sec,  depending  on 
the  desired  resolution.  Thus,  the  results  of  this  portion  of  the  analysis  may  be 
directly  compared  with  other  time-based  regular  data  bases  such  as  the  Unified 
Abstract  (UA)  files  from  the  Atmosphere  Explorer  project.  The  results  of  this 
analysis  are  especially  useful  in  the  very  low  attenuation  region  where  the  standard 
data  production  program  LMC80  is  not  particularly  useful.  For  these  cases  the 
data-fitting  range  of  LMC85  may  be  expanded  and  the  form  of  the  linear  regression 
for  low  attenuation  usually  matches  the  true  change  in  the  date  more  closely  than 
does  the  exponential  regtession  used  by  LMC80.  In  addition.  LMC85  incorporates 
model-based  corrections  for  data  on  Sq  taken  from  a  region  where  the  signals  are 
significantly  attenuated,  in  order  to  provide  a  better  value  of  optical  depths  cal¬ 
culated  with  the  so  corrected  value  of  the  unattenuated  reference  signal  SQ.  Another 
important  application  of  the  program  LMC85  is  its  occasional  use  for  diagnostic 
comparison  of  its  results  with  those  obtained  by  the  standard  program  L.MC80.  This 
provides  a  valuable  check  on  the  results  of  both  systems,  as  they  both  operate  on  the 
same  observational  data  base. 
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4.2  IMmnuution  ut  "l  mttcniulrd'  Hrlmnrr  Sjjiul  (S(>) 

4.  2.  1  KANGE  SELECTION 

As  in  the  standard  data  reduction  system  (LMC80),  the  reference  signal  for  each 
wavelength  processed  by  the  special  program  LMC85  may  be  determined  in  several 
ways;  however,  only  one  S^-eligible  region  is  used  for  each  wavelength,  and  all  of 
the  wavelengths  using  a  specific  method  also  share  a  common  S^-eligible  region. 
Each  wavelength  may  have  its  reference  signal  computed  by;  (a)  direct  entry  of  the 
raw  reference  signal,  (b )  automatic  selection  of  an  SQ-eligible  region,  or  (c)  entry 
of  the  time  land  time  span)  of  an  "S^-eligible  region"  from  which  the  reference 
signal  will  be  extracted.  For  the  entry  of  an  S^-eligible  region,  the  experimenter 
enters  the  date  and  time  of  the  center  of  the  chosen  region  as  well  as  the  size  of  the 
time  interval  for  the  corresponding  data  fit.  The  resulting  processing  time  interval 
then  includes  the  given  interval  before  the  given  time  and  the  interval  after  it.  All 
wavelengths  for  which  this  method  has  been  chosen  use  this  same  interval.  Typical 
intervals  over  which  S  is  computed  are  of  the  order  of  15  sec,  corresponding  to  a 
total  SQ-eligible  region  of  30  seconds.  The  interval  used  is  carried  into  the  file 
along  with  an  indication  of  the  method  of  S  - range  selected.  For  the  direct  entry 
of  a  raw  unattenuated  signal,  a  date  and  time  must  be  included  for  use  in  obtaining 
the  orbit  data  for  the  absorption  correction  described  in  Section  4.  2.  3.  Each  wave¬ 
length  with  the  direct  entry  of  a  raw  unattenuated  signal  can  have  a  different  date 
and  time.  The  final  method  of  S^-eligible  region  selection  (automatic  selection  by 
the  program)  is  accomplished  by  selecting  the  time  of  minimum  solar  zenith  angle. 
Since  this  section  of  the  system  is  used  primarily  for  near  equatorial  circular  orbit 
conditions,  the  region  of  least  attenuation  can  be  simply  selected  by  this  method. 
There  is  no  attempt  to  ensure  a  maximum  tolerance  of  attenuation,  since  all  S 

o 

values  receive  a  correction  for  this  (see  Section  4.  2.  3)  and  the  magnitude  of  this 
correction  is  carried  into  the  resulting  data  files.  The  exact  SQ-time  for  the  fitted 
value  is  then  taken  as  the  nearest  15-sec  point  to  the  time  of  minimum  solar  zenith 
angle.  The  size  of  the  time  interval  is  determined  the  same  way  as  for  a  direct 
entry.  Again,  the  given  time  range  (in  seconds)  is  taken  both  before  and  after  the 
center  of  the  region.  If  the  center  point  of  the  region  is  less  than  the  given  interval 
from  either  end  of  the  pass,  or  if  the  EUVS  instrument  is  scanning  for  a  flux  study 
at  that  time,  the  S  -eligible  region  is  repositioned  accordingly,  again  to  a  whole 
15-sec  point.  As  is  the  case  for  the  direct  entry  method,  all  wavelengths  for  which 
this  method  has  been  specified  will  share  the  same  S^-eligible  region,  but  this 
region  need  not  be  the  same  as  that  of  the  wavelengths  processed  with  the  direct - 
entry  method,  nor  any  >>f  those  wavelengths  where  the  raw  signal  and  a  reference 
date  and  time  have  bee'  given  by  the  experimenter. 
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4.2.2  KEG  KESSION  ALGORITHM 


The  value  of  the  raw  unatteauated  signal  for  each  wavelength  is  determined  by 
a  standard  linear  least  squares  regression  of  the  form 

S  =  A  +  B  At  .  (9) 


In  this  equation,  S  is  the  raw  signal,  At  is  the  elapsed  time  in  seconds  from  the 
start  of  the  fitting  interval,  and  A  and  B  are  the  regression  coefficients.  The  same 
checks  for  the  minimum  number  of  points  and  for  the  insignificant  determinants  are 
performed  as  in  the  standard  data  reduction  system  (LMC80).  The  unattenuated 
raw  signal  is  taken  from  the  center  of  the  fitting  interval,  and  the  standard  error 
is  computed  and  kept  in  the  same  record  with  the  data. 

4,2.3  FINAL  S0  DETERMINATION  AND  BACKGROUND 
COR  RECTION 

Once  the  raw  unattenuated  signal  is  extracted  for  each  wavelength,  it  is  cor¬ 
rected  for  background  bv  the  same  algorithm  as  was  used  in  the  regular  data  pro¬ 
duction  system  (see  Section  3.2.3).  However,  after  this  has  been  done,  the  result¬ 
ing  signal  is  corrected  for  residual  attenuation  by  the  atmosphere.  This  correction 
is  obtained  by  the  use  of  one  of  the  models  mentioned  below  to  predict  particle 
number  densities  of  O,  O0.  N.,,  and  He.  From  these,  constituent -associated  optical 
depths  are  computed  and  summed  to  find  a  predicted  total  value  of  optical  depth  for 
"SQ.  "  Once  this  optical  depth,  r  has  been  obtained,  the  corrected  value  of 

(S  )  .  is  obtained  from  the  relation 

o  corr 


(S  )  =  S  e 

o  corr  o 


(10) 


Three  different  atmospheric  models  may  be  used  to  generate  the  predicted  value  of 

3  4 

T  (S  ).  that  is,  the  MSIS  model,  see  Hedin  et  al.  ’  the  Jacchia  1971  model,  or  a 
o  o 

default  model  which  uses  an  input  temperature  and  set  of  densities  at  a  given  refer¬ 
ence  height  along  with  an  assumed  exosphere  temperature  and  "shape  factor"  to 
generate  the  final  densities  using  a  "Bates  profile.  "  The  constituent  optical  depths 
and  an  indication  of  the  model  used  to  compute  them  are  carried  into  the  output,  so 

3.  Hedin,  A.  K. ,  Salah.  J.  E.  ,  Evans.  J.  V.  .  Reber,  U.  A.  ,  Newton,  G. 

Spencer.  N.  W. ,  Kayser,  D.  C. ,  Alcayde,  C. ,  Bauer,  P. ,  Cogger,  L. ,  and 
McClure,  J.P.  (1977)  A  global  thermospheric  model  based  on  mass  spectrom¬ 
eter  and  incoherent  scatter  data,  MS1S  1.  N„  density  and  temperature. 
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that  the  correction  may  be  removed  and  another  one  substituted  using  a  different 
model  whenever  this  may  become  desirable  at  some  later  time.  For  most  pur¬ 
poses.  the  MS1S  model  was  used,  and  the  default  model  was  used  only  for  testing 
purposes.  Although  the  model  used  has  errors  in  the  predicted  number  densities 
that  can  be  quite  large,  the  effect  upon  the  i>o  data  is  usually  small.  For  S  regio 
which  are  nearly  unattenuated,  the  predicted  optical  depths  will  be  so  small  that  t 


In  this  equation,  n.  is  the  predicted  particle  number  density,  o.  is  the 

cross  section  for  the  wavelength  being  corrected  for  that  constituent 

partial  pressure  scale  height  of  that  constituent,  and  F.  is  the  so  cal 

depth  factor  required  to  convert  the  vertical  optical  depth,  T  (T 

v,]  v,. 

the  optical  depth  T  ^  for  the  given  slant  column,  (T  .  =  Ty  .  F  ).  The 
11.  are  computed  from  the  relation 


(15) 


For  x -  105° 


,  In  He+hq 

j  =  V5  MJ 


,  H. 

°  +  8  Ri) 
e 


1  h  -  h. 

1  4  erf  J 

s  Q 

ii. 
j 

These  equations  assume  a  spherical  atmosphere,  and  compute  a  correction  to  the 
optical  depth  at  the  point  of  minimum  ray  height,  h„.  The  point  Q  is  called  remote 
probing  point  for  x  >  90°,  whereas  Q  is  identical  with  the  actual  satellite  point  for 
X  -  90°  (see  Figure  1), 

Upon  the  completion  of  the  above  correction  process,  (summing  T  .  for  all 
constituents),  the  value,  its  date  and  time,  its  corresponding  orbit  data,  and  the 
experimenter  input  controls  are  written  to  the  output  file.  All  of  the  orbit  data  used 
in  the  computations  and  for  output  have  been  corrected  for  possible  instrumental 
misalignment  by  a  process  described  in  detail  in  Section  4.5. 


4.5  Opliral  Depth  ('.ompiilatioiii* 

4.  3.  1  RANGE  SELECTION 

Optical  depth  data  are  computed,  regressions  performed,  and  the  resulting 
fitted  values  at  the  requested  time  intervals  are  reported,  along  with  associated 
orbit  parameters.  The  format  of  each  entry  generated  in  the  data  file  is  specified 
in  Table  5.  There  can  be  any  number  of  these  entries  following  the  entry  for  S 
and  its  associated  parameters.  Two  time  intervals  in  sec  are  given  bv  the  user  of 
the  program  LMC85.  The  first  is  the  data-taking  interval,  T.  The  system  takes 
data  starting  from  the  given  starting  date  and  time  through  the  ending  date  and  time 
once  in  each  interval.  These  times  are  synchronized  to  0  sec  on  the  starting  date, 
that  is,  the  first  point  occurs  for  a  time  t  =  N.  T,  where  N  is  the  lowest  integer 
value  that  will  give  data  within  the  pass.  The  volume  of  data  received  depends  solely 
upon  this  parameter.  Intervals  which  are  fractions  or  multiples  of  15  sec  are  par¬ 
ticularly  useful,  as  they  enable  direct  comparisons  of  EUVS  data  with  the  Atmosphere 
Explorer  Unified  Abstract  data  base.  The  actual  parameter  used  for  any  given 
reduction  is  given  in  the  last  section  of  the  S  entry  for  that  pass.  This  interval  is 
the  same  for  all  wavelengths  under  analysis  and  has  no  effect  upon  the  length  of  the 
fitting  interval  which  may  overlap  the  adjacent  ones.  The  latter  is  determined  by 
the  second  time  interval  given  by  the  program  user  and  it  works  in  the  same  way  as 
that  for  the  S  -  computation.  This  means  that  all  data  from  that  interval  before,  to 
the  interval  after,  are  included  in  the  regression.  The  latter  interval  selection  is 
also  reported  in  the  S(j-associatcd  entry.  It  is  independent  of  the  corresponding 
interval  selected  for  the  Sq- computation. 
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fable  3.  File  I'oniml  for  Optical  Depth  Special  Studies  Files 


Flag  indicating  nod  4  records  of  the  format,  and  nil 

subsequent  optical  depth  records  will  refer  to  if. 

S  Krooni  1 

S  Kesults 

o 

O 

i'«.  a 

S  (1'ouutsl 

f«.  a 

llackg  round  (Founts! 

F5.  .1 

Optical  depth  correction  for  constituent  O 

F5.  3 

Optical  depth  correction  for  constituent 

F9.  3 

Optical  depth  correction  for  constituent  Oy 

F5.  3 

Optical  depth  correction  for  constituent  Me 

A  1 

Satellite  in  (C,  1',  or  K) 

19 

Mate  of  S  (Julian  da\  vvddd) 

Id 

l  ime  of  S  (milliseconds! 

A  4 

Method  of  S  (  eligible  region  range  selection 

Fl’Tl  automatic  (lowest  solar  zenith  angle! 

KN  1  It  -  entry  of  region 

DIPT  -  direct  entry  of  raw  value  and  time 

S  Hfcorii  2 
o 

S  IfrAulta 

0 

1 PK 10.  3 

.'eiv  order  regression  coefficient  (  VI 

ll’K  10.  3 

First  onler  regression  coefficient  (10 

OPF6,  a 

Standant  Krmr  of  concession 

13 

Number  of  points  in  res  cession 

13 

Number  of  upper  limit  violations 

13 

Number  of  lower  limit  violations 

S  Rt>eord  3 

-  Orbit  Parameters 

19 

Orbit  Number 

re.  a 

Solar  Zenith  Annie  (den' 

F7.  a 

Minimum  Kay  Height  (km! 

F7.  a 

Satellite  Height  (km) 

Ft.  a 

Kemote  Probing  Point  1  atltudc  (deg! 

i'7.a 

Kemote  Probing  Point  Fast  1  ongitude  (deg! 

F6.  a 

Kemote  Probing  Point  Mlpole  1  atttude  (deg! 

Id 

Kemote  Probing  Point  1  oral  Solar  I'tme 
(hlunmss diours  minutes  seconds' 

u; 


t 


Table  5.  File  Format  for  Optical  Depth  Special  Studies  Files  (Conti 


S  Record  4 
o 

-  F.  q  Ap,  and  input  parameters 

F5.  1 

Fj0  ,  Daily  Average  for  Sq  date 

F5.  1 

Fj()  _  Daily  Average  for  Day  prior  to  SQ  date 

13 

Average  Ap  value  for  9  to  6  hr  prior  to  S^  time 

13 

Average  AP  value  for  6  to  3  hr  prior  to  Sq  time 

13 

Average  Ap  value  for  3  to  0  hr  prior  to  S  time 

13 

Average  Ap  value  for  3  to  0  hr  prior  to  S  time 

16 

Reference  file  key  ID  (yvmmdd) 

16 

Model  and  Cross  Section  Reference  file  key  ID  (yvmmdd) 

16 

Date  of  System  Execution  (yvmmdd) 

F4.  1 

Range  (in  sec)  around  S  center  data  and  time 

F4.  1 

Range  (in  sec)  around  optical  depth  date  and  time 

F4.  1 

Data  taking  interval  for  optical  depths  (sec) 

Each  Optical 

Depth  Entry  comprises  two  records  of  the  following  format* 

and  are  indefinite  in  number  for  each  S  entry. 

o 

Optical  Depth  Record  1  -  Optical  Depth  and  Regression  Results 

15 

Date  of  Optical  Depth  Point  (yyddd) 

18 

Time  of  Optical  Depth  Point  (msec) 

1*’7 .  3 

Optical  Depth  Value 

1PE10.  3 

Zero  order  regression  coefficient 

1PE10.  3 

First  order  regression  coefficient 

OPF6.  2 

Standard  error  of  regression 

13 

Number  of  points  in  regression 

13 

Number  of  upper  limit  violations 

13 

Number  of  lower  limit  violations 

Optical  Depth  Record  2  -  Optical  Depth  Orbit  Parameters 

F6.  2 

Solar  Zenith  Angle  (deg) 

F7.2 

Minimum  Rav  Height  (km) 

F7.  2 

Satellite  Height  (km) 

F6.  2 

Remote  Probing  Point  Latitude  (deg) 

F7.  2 

Remote  Probing  Point  East  Longitude  (deg) 

F6.  2 

Remote  Probing  Point  Dipole  Latitude  (deg) 

16 

Remote  Probing  Point  Local  Solar  Time 

(hhmmssthours  minutes  seconds) 
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4.8.2  REGRESSION  ALGORITHM 


Optical  dcptli  data  are  obtained  from  standard  linear  least  squares  regressions 
of  the  form 

T  =  A  +  B  At  ,  (16) 

In  this  equation,  At  is  the  time  in  elapsed  seconds  since  the  stai-t  of  the  fitting 
interval,  and  A  and  B  are  regression  coefficients.  As  in  the  standard  program 
(LMC80),  the  optical  depths,  T,  at  each  point  are  determined  as 

r  =  -  Cn <SQ/S)  ,  (17) 

where  Sq  is  the  corrected  unattenuated  signal  discussed  above.  The  same  checks 
for  a  minimum  number  of  data  points  or  singular  system  are  made  as  in  the  S^ 
regression,  and  the  standard  error  is  also  computed  and  reported  in  each  optical 
depth  entry  made  in  the  resulting  data  files. 

4.4  Data  filtering 

In  keeping  with  the  special  nature  of  the  program  LMC85,  the  upper  and  lower 
limits  for  data  acceptance  are  left  somewhat  wider  than  they  are  for  the  standard 
program  LMC80.  Only  for  the  Sq  computations,  the  data  filtering  is  performed  in 
the  same  way  as  in  LMC80,  that  is,  reference  values  of  upper  and  lower  limits  are 
used  for  each  wavelength.  Afterwards,  in  the  optical  depth  computations,  the  lower 
limit  is  set  to  zero,  and  the  upper  limit  is  set  to  a  value  given  as 

S  ..  -  S  +  5  Js~  .  (in) 

upper  limit  o  }/  o  'io' 

This  method  is  used  to  eliminate,  not  merely  so-called  fill  data,  but  also  data  which 
are  grossly  aberrant,  for  example,  due  to  transmission  errors  of  "stuck-bit" 
problems. 

As  was  the  case  in  the  standard  data  production  (LMC80),  the  magnitude  of  the 
instrumental  count  sample  for  S  determines  the  practical  maximum  value  of  optical 
depths  that  can  be  evaluated.  Again,  data  for  optical  depths  corresponding  to  a 
signal  of  less  than  five  counts  (after  background  correction)  will  not  be  placed  in 
the  output  files. 
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4.5  Alignment  Correction*  to  Orbit  l>*ta 

4.5.1  HATIONALE 

The  EUVS  instrument  is  pointed  by  the  Solar  Pointing  System  (SPS)  at  any  one 
of  256  different  target  points  covering  a  square  angular  field  slightly  greater  than 
the  solar  disk.  For  routine  geophysical  observations,  the  target  is  selected  close 
to  the  center  of  the  solar  disk.  Due  to  finite  mechanical  tolerances,  there  are  some 
slight  departures  from  exact  co-alignment  for  each  of  the  24  monochromators  of 
the  EUVS  instrument.  As  a  result,  each  monochromator  entrance  axis  has  some 
slight  angular  departure  (on  the  order  of  several  minutes  of  arc)  from  the  axis  of 
SPS  pointing,  and  thus,  for  narrow-field  monochromators,  a  slightly  different 
effective  value  of  solar  zenith  angle.  These  small  differences  in  the  solar  zenith 
angle  lead  to  correspondingly  small  changes  in  the  minimum  ray  height  (h.  )  at  which 

•  W 

data  are  taken  by  each  monochromator  (see  Figure  1).  This  effect  is  negligible  for 
those  monochromators  with  full  disk  view  of  the  sun.  Even  for  narrow-field 

monochromators,  these  departures  are  quite  significant  for  large  values  of  h.. 

W 

(where  the  scale  height  is  large).  However,  for  high  solar  zenith  angles,  the  effec¬ 
tive  value  of  the  remote-probing  height  h^  from  that  of  perfect  coalignment  may  be 
several  kilometers.  This  effect  is  obviously  more  serious  for  wavelengths  pene¬ 
trating  deeply  into  the  atmosphere,  for  example,  H  Lv-a  ,  where  the  scale  heights 
are  low  and  the  absorption  increases  rapidly  once  attenuation  begins.  This  problem 
is  handled  by  computing  corrected  solar  zenith  angles  and  minimum  ray  heights 
separately  for  each  monochromator,  that  is,  for  each  wavelength.  All  other  orbit 
parameters  included  in  the  data  file  (latitudes,  longitudes,  local  solar  time,  and 
so  on)  are  also  computed  for  these  corrected  values  of  solar  zenith  angle,  although 
the  small  differences  in  the  latter  are  usually  not  significant. 

4.  5.  2  ALGOKITHM 

In  order  to  compute  the  true  solar  zenith  angle  for  those  wavelengths  which  can 
have  a  significant  deviation,  a  data  set  of  offsets  is  provided  (one  each  for  the  t  - 
and  y-  directions  respectively)  for  each.  These  offsets  were  determined  from 
occasional  diagnostic  modes  of  operation  called  "raster  scan,  "  where  the  SPS  target 
point  rasters  across  the  solar  disk.  An  evaluation  of  the  responses  of  the  various 
monochromators  during  such  a  raster  then  allows  one  to  determine  the  offsets  of 
each  monochromator-entrance  axis  with  respect  to  the  ideal  point  direction.  The 
resulting  offsets  expressed  in  SPS  steps  were  then  converted  to  angular  offsets  in 
the  €  -  and  y  -  directions  and  placed  in  a  table  for  use  by  the  program  LMC85  and 
other  programs  (see  Table  6). 

Detailed  discussion  on  the  correction  for  solar  alignment  of  monochromators 
and  the  effects  of  fields-of-view  are  shown  in  Appendix  A,  Sections  A1  and  A2, 
respectively. 


Table  6.  Alignment  Offsets  for  AE-C  in  Kadians 


Wavelength 

Epsilon 

Gamma 

1 

175  A 

2 

465 

3 

584 

No  correction 

required 

4 

765 

5 

855 

6 

1609 

-0. 00206 

+0. 00124 

7 

1457 

-0,  00206 

+  0.  00124 

8 

1700 

-0,  00206 

+  0.  00124 

9 

1730 

-0, 00206 

+  0,  00124 

10 

1835 

-0.  00206 

■+  0.  00124 

11 

256 

-0, 00204 

+0. 00162 

12 

304 

-0. 00204 

+  0.  00019 

13 

610 

-0.  00204 

+  0.  00117 

14 

465 

-0. 00204 

+  0.  00099 

15 

584 

-0.  00219 

+0. 00158 

16 

1026 

-0. 00198 

-0. 00147 

17 

977 

-0.  00204 

-0. 00195 

18 

1216 

-0.  00097 

-0, 00007 

19 

1600 

-0. 00369 

+0. 00022 

20 

1775 

-0. 00369 

+0. 00047 

NOTE:  For  AE-C  orbits  after  860,  the  EUVS  instrument 
was  offset  7  arc  min  in  the  positive  epsilon 
direction  to  improve  the  pointing  of  the  narrow 
field-of-view  monochromators,  that  is  to  have 
them  point  more  closely  to  the  center  of  the 
solar  disk.  Therefore,  this  offset  (0.  0020362175 
radians)  was  added  to  the  epsilon  offset  for 
these  orbits. 

5.  MODEL  COMPARISON  (PROGRAMS  LMCURSI  and  LMCIAGA) 

5.1  Introduction 

The  purpose  of  the  model-comparison  programs  (LMCUltSI  which  accesses 
LMC80  produced  data  and  LMCIAGA  which  accesses  LMC85  produced  data)  is  to 
provide  a  comparison  of  several  atmospheric  models  with  the  EUVS-derived  results 
on  optical  depths  for  varying  conditions.  To  accomplish  this,  the  model-predicted 
particle  number  densities  for  atomic  oxygen,  molecular  nitrogen,  molecular  oxygen 
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and  helium  are  converted  to  both  a  total  optical  depth  and  a  three-constituent  optical 
depth  (adding  contributions  by  O,  N£.  and  Og  only)  using  a  standard  data  set  of 
absorption  cross  sections.  The  ratios  of  these  model-values  of  optical  depths  to 
the  values  derived  from  the  EUVS  observations  are  then  determined.  These  ratios, 
along  with  all  other  associated  details  including  identification  of  model  and  list  of 
all  relevant  geophysical  parameters  for  each  point  in  time,  are  then  sent  to  output 
files  in  the  EBCDIC  format.  In  addition  to  the  ratios,  each  entry  contains  the  date 
and  time,  the  constituent  optical  depths  generated  from  the  selected  model,  the 
total  optical  depth  from  EUVS,  the  helium  contributed  optical  depth  used  with  it, 
and  the  associated  orbit  parameters.  An  example  of  the  actual  output  of  the  pro¬ 
gram  LMCURSI  is  given  in  Table  7. 

Any  one  of  three  models  may  be  selected  by  the  experimenter,  and  the  result- 
ing  file  entry  indicates  which  model  was  selected.  These  are:  the  MS1S  model  ; 
the  1971  Jacchia  model  (C1KA  1972);  or  a  so  called  "default"-model  which  is  a 
simple  test  model  iti  which  the  user  of  the  programs  selects  a  reference  altitude, 
the  constituent  particle  number  densities,  and  temperature  at  that  reference  alti¬ 
tude  as  well  as  a  shape  factor  (using  a  "Bates"  profile)  and  an  exosphere  tempera¬ 
ture.  These  input  parameters  are  then  used  by  the  default  model  to  generate  the 
constituent  densities  and  scale  heights  at  any  other  altitudes  for  which  model-based 
values  of  optical  depths  are  desired.  It  should  be  noted  that  this  default  model 
provides  no  adjustments  for  diurnal,  latitudinal,  seasonal,  or  other  types  of  atinos- 
spheric  variabilities.  This  means  that  this  default  model  is  strictly  an  ad  hoc 
utility  allowing  experimentation  with  various  sets  of  model  inputs  such  as  may  be 
desired  for  a  detailed  study  of  one  specific  set  of  observations. 

.1.2  Optical  Depth  Ecru-ration  ami  Model  Oonskleralions 

Although  the  programs  LMCURSI  and  LMCIAGA  produce  a  similar  result,  they 
operate  on  different  data  bases. 

LMCURSI  uses  so-called  "GU  files,"  that  is,  the  binary  files  generated  by  the 
standard  program  LMC80.  Due  to  the  compact  form  required  for  these  GU  file 
records,  they  have  none  of  the  associated  orbit  parameters  and  they  do  not  provide 
regular  coverage  in  time.  However,  due  to  the  compactness  of  those  GU  files  and 
the  regular  coverage  provided,  a  wide  range  of  geophysical  conditions  can  bo 
explored  in  a  single  execution  of  LMCURSI. 


Since  many  of  the  SQ-eligible  regions  were  directly  input,  that  is,  selected  in 

spite  of  significant  residual  optical  depth  to  increase  the  rate  of  data  return,  a 

correction  for  these  residual  optical  depths,  To  (Sq),  is  performed  at  this  level. 

The  optical  depth  of  each  date  and  time  is  computed  from  the  model-predicted 

densities  and  temperatures  by  using  the  same  algorithm  described  in  Section  4.  2.  3 

and  is  used  as  the  T  (S  ).  That  is,  the  value  of  T  (S  )  is  added  to  the  standard 
o  o  o  o 

optical  depth  values  of  each  GU  file  entry  within  that  pass.  Thus  the  program 
LMCUHSI  generates  values  of  corrected  optical  depths  corresponding  to  more  truly 
unattenuated  SQ  values.  For  double  passes,  where  two  S  values  are  available, 
both  residual  optical  depths  are  computed,  and  the  final  r  ^  used  for  each  optical 
depth  point  is  the  linearly  interpolated  value  between  them.  For  the  sake  of  con¬ 
sistency,  the  same  model  is  used  for  the  T  ^  computation  as  that  requested  for  the 
ratio  comparison. 

Since  the  LMC85-created  special  study  files  already  contain  a  residual  optical 
depth  correction,  the  model-comparison  program  LMC1AGA  operating  on  these 
files  first  removes  the  LMC85-inserted  model  correction,  and  in  its  place,  a  new 
correction  is  performed  with  that  model  which  is  specified  for  the  present  execution 
of  LMC1AGA. 
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Appendix  A 


Correction  for  Solar  Alignment  of  Monochromators 
and  Effects  of  Fields-of-View 


Al.  CORRECTIONS  FOR  SOLAR  ALIGNMENT  OF  MONOCHROMATORS 


The  solar  pointing  subsystem  (SPS)  establishes  a  reference-target  point  on  the 
solar  disk  with  an  accuracy  of  about  0.0003  radian.  Figure  Al  gives  a  schematic 
overview  of  geometrical  relationships.  The  actual  dimensions  of  the  fields -of -view 
of  the  so-called  "narrow -field  monochromators"  are  still  relatively  large,  for 
example,  0.001  radian  for  the  monochromators  No.  12,  No.  23,  and  No.  24  in  both 
elevation  (f  )  and  cross -elevation  (y)  and  range  up  to  0.003  radian  lor  some  other 
narrow -field  MN's  of  the  EUVS  instrument.  Detailed  illustrations  of  all  relevant 
geometrical  aspects  are  given  in  Figures  A2,  A3,  and  A4.  Therefore  the  absolute 
SPS-associated  pointing  accuracy  is  indeed  quite  adequate.  However,  some  of  the 
narrow -field  MN's  on  the  AE-C  satellite  were  found  to  have  individual  entrance  axes 
departing  from  the  SPS-pointing  direction  by  angles  ranging  up  to  0.0036  radian  in 
£  direction  and  about  0.  002  radian  in  the  y  direction  (see  Figures  Al  and  A2)  for 
the  worst  cases.  The  actual  values  of  the  MN-peculiar  misalignment  angles  have 
been  determined  by  evaluation  of  data  obtained  in  diagnostic  modes  of  operation  of 
the  SPS,  that  is,  from  16-step  e  -scans,  16-step  y -scans  and  256-step  full  raster 
scans.  The  accuracy  of  this  determination  is  estimated  to  be  equal  to  that  of  the  SPS 
accuracy,  that  is,  also  about  0.  0003  radian.  Within  that  accuracy,  no  changes  of 
these  misalignment  angles  were  expected  (considering  the  design)  and  none  were 
indeed  observed.  The  accuracies  quoted  had  to  be  accepted  as  adequate  for  the 
objectives  of  aeronomical  EUV  absorption  analysis  planned  for  the  EUVS  experiment. 
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even  though  they  would  of  course  be  considered  unacceptably  poor  for  instruments 
with  truly  narrow  fields-of-view. 

Whereas  the  selection  of  fields-of-view  of  less  than  about  1  arc  min  would  have 
been  clearly  undesirable  for  aeronomical  EUV  absorption  analysis  (see  Garriott 
et  al)'^,  fields  significantly  narrower  than  those  of  the  EUVS  instrument  for  AE 
would  of  course  have  been  desirable,  along  with  a  better  accuracy  of  monochromator 
co-alignment  and  somewhat  better  absolute  accuracy  of  verifying  SPS- pointing 
directions.  However,  the  actual  EUVS-design  choices  had  to  be  made  within  rather 
strict  constraints  of  cost,  size,  weight,  power  consumption,  and  lead-time  to  com¬ 
pletion  of  instruments  for  flight,  aside  from  the  emphasis  on  full-disk  measure¬ 
ments.  These  constraints  necessitated  certain  otherwise  undesirable  compromises 
the  consequences  of  which  are  analyzed  in  Sections  A1  and  A2.  The  latter  deals 
with  the  problems  due  to  the  relatively  large  fields-of-view  per  se,  that  is,  prob¬ 
lems  which  would  remain  even  if  co-alignment  and  SPS  pointing  had  been  absolutely 
perfect. 

The  situation  of  monochromator  co-alignment  in  the  EUVS  instruments  on  AE-D 
(the  spacecraft  failed  in  January  1976)  and  AE-E  (still  operating)  was  significantly 
better  than  that  for  AE-C  (operated  until  December  1978).  Only  for  the  AE-C 
EUVS,  the  departures  of  MN-peculiar  pointing  directions  from  the  SPS-established 
reference  direction  exceed  the  possible  error  of  experimental  verification  sig¬ 
nificantly.  The  AE-C  monochromators  also  showed  a  significantly  asymmetric 
pattern  of  target  points  around  the  SPS-established  reference  target.  Therefore, 
geophysical  routine  observations  from  AE-C,  originally  planned  to  be  made  with 
center-pointing  of  the  SPS  (NE  =  NG  =  8  in  Figure  Al),  were  changed  to  use  a 
standard  offset  of  two  elevation  steps,  that  is,  NE  =  10,  NG  =  8  in  Figure  Al,  re¬ 
sulting  in  values  of  t  <,p<.  =  0.  002036  radian  and  =  0.  For  this  mode  of  opera¬ 

tion,  retained  for  all  geophysical  observations  after  AE-C  orbit  860  (March  1974), 
the  pattern  of  MN-peculiar  target  points  characteristic  of  the  AE-C  instrument 
represents  optimum  symmetry  around  the  center  of  the  solar  disk  and  also  shows 
significantly  reduced  extreme-values  of  MN-pointing  departures  from  the  solar 
vector,  that  is,  none  in  excess  of  about  0.0015  radian.  For  both  AE-D  and  AE-E, 
the  corresponding  optimum  was  indeed  that  of  no  offset  from  sun-center  pointing  of 
the  SPS,  with  individual  MN-misalignment  angles  eventually  no  greater  than  the 
estimated  accuracy  of  their  verification.  Therefore,  Section  Al  is  important  for 
AE-D  and  AE-E  data  only  because  of  its  use  in  the  evaluation  of  effects  of  the  finite 
fields-of-view— discussed  in  Section  A2. 

Our  correction  of  effective  solar  zenith  angles  and  remote -probing  heights, 
accounting  for  the  known  departures  of  the  individual  monochromator-entrance  axes 


5.  Garriott,  O.K.,  Norton,  R.  B. ,  and  Timothy,  J.G.  ( 1977)  Molecular  oxygen 
concentration  and  absorption  cross  sections  in  the  thermosphere  derived 
from  extreme  ultraviolet  occultation  profiles,  J.  Geophys.  Res. 

82(No.  32):4973-4982. 

ww 


38 


i 


from  .he  direction  of  the  «l„  vector  I?  t.  detined  to  (rom 

position  the  center  „  ,h.  eoi.r  dish,  .'t,  hosed  on  th.  r.l.tZ 

y,  -  (s,  X  a,  J  /  |  *I  x  af  | 


1  I  ‘  lai  '  si  (si  *  »!>]  /  Isj  X  Tjl 


<Al) 


where  the  geometry  of  these  unit  vectors  e  y  s  and  *  „ 
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The  actual  amount  of  data  rejected  by  the  exclusion  of  those  for  G  i  87*  has 
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a  j  =  n}  (aj  •  s  j )  /  I  tij  •  Sj|  =  «j  sign  (n,  •  Sj) 

may  appear  crude,  but  the  actual  consequences  for  our  purposes  of  calculating  e  } 
and  y  j  from  the  relations  [Eq.  (1))  are  readily  acceptable.  The  angles  formed  by 
the  directions  of  t*j  and  y  t  with  some  fixed  reference  direction  in  the  plane 
perpendicular  to  the  solar  vector,  for  example,  the  direction  of  the  projected  solar 
axis  of  rotation,  do  not  need  to  be  known  very  accurately  — say  to  a  few  degrees. 

For  instance,  a  rotation  of  ?  and  y  j  by  a  full  5°  in  this  plane  which  is  the  plane  of 
drawing  in  Figure  Al(a),  would  change  the  offset  coordinates  of  a  nominal  target 
point  only  by  about  0.  0001  or  0.  0002  radian  if  that  target  point  is  offset  from  null 
by  0.  001  or  0.  002  radian  respectively,  corresponding  to  1  or  2  steps  of  SPS-offset. 
Since  none  of  the  MN-misalignment  angles  relative  to  the  solar  vector  were  greater 
than  0.  003  radian  even  for  the  worst  case  of  the  AE-C  MN's,  the  adoption  of  rela¬ 
tion  [Eq.  (A2)J  is  indeed  justified  as  long  as  one  excludes  the  region  of  angles 
G  <  7*.  for  which  the  instrumental  accuracy  of  offset  pointing  by  the  SPS  also 
begins  to  become  marginal. 

Defining  the  angle  E  as  the  elevation  of  the  sun  measured  from  the  satellite's 
orbit  plane  and  defining  the  sign  of  E  as  positive  (negative)  when  the  sun  is  on  the 
positive  (negative)  side  of  the  orbit  plane,  one  finds  sin  E  =  n j  •  Sj  and  can  write 
the  relations  [Eq.  (AD)  with  adoption  of  the  approximation  [Eq.  (A2)]  as 


y  =  ±  Sj  X  Sj/cos  E  , 
t  j  =  ±  n  j  /  cos  E  +  s  j  |  tan  E  |  , 


(A3) 


where  the  upper  or  lower  signs  must  be  used  for  E  >  0  or  E  >  0,  respectively. 

Defining  the  instrumental  target  point  on  the  solar  disk  as  the  center  of  the  area 
of  the  solar  disk  viewed  by  a  given  monochromator,  one  finds  the  unit  vector  of  the 
direction  from  the  satellite  to  that  target  point  (briefly  called  "target  vector"  below) 
to  be  given  as 

=  (Ij  +  e£j  +  r?i)  (i  +  t2  +  y2)  ^2.  <A4) 


where  t  =  t '  +  tgpg  and  y  =  y  '  +  I'gpg  are  the  offset  coordinates  of  the  instrumental 
point  of  a  given  monochromator  relative  to  the  center  of  the  solar  disk;  t  '  and  y  '  are 
the  corresponding  offsets  relative  to  the  SPS-established  reference  target;  tgpg  and 
are  the  coordinates  of  the  SPS  target  relative  to  the  center  of  the  solar  disk. 
This  expression  of  the  target  vector  [Eq.  (A4)]  in  connection  with  the  relations 
[Eq.  (A3)]  defining  the  chosen  offset-coordinate  axes,  has  been  used  in  all  of  our 
computer  programs  where  solar  zenith  angles  (x)  and  re  mote -probing  heights  (h^) 
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are  required  to  bo  referred  to  that  particular  .solar  ray  which  connects  the  center 
of  the  instrumentully  viewed  area  of  the  -solar  disk  with  the  point  of  observation 
(satellite).  It  should  bo  noted  that  the  same  relations,  jEq.  (Ad))  and  | lSq.  (A4)|, 
also  serve  as  tin-  basis  for  the  assessment  of  the  effects  of  finite  fields  of  view 
discussed  In  Section  A2, 

The  zenith  angle  associated  with  a  given  target  point  is  then  determined  by  the 
relation 

-%  a 

COS  x-  P,  •  tt  (A5) 

where  |V  ■  1^/1*  (with  1*  *  |1^|  )  is  the  unit  vector  of  the  direction  from  the  earth's 

1  4 

center  to  the  satellite  and  t(  is  the  target  vector  (also  unit  vector)  introduced  above. 

The  "remote-probing  point"  associated  with  that  target  point  on  the  solar  disk 
is  then  determined  by  the  geocentric  position  vector  given  as 

Q  -  P  -  tjp  cos  x;  Q  (\  S  ilO")  ■  p  ,  (AG) 

and  the  geocentric  distance  of  the  probing  point  is  given  by  its  absolute  magnitude, 

Q  *  |Q|  ■  P  sin  x:  W(X  1  90*)  «  P  .  (A7) 

The  "probing  height"  is  obtained  as 

"  «  *  »E<0g>i  hg  <X  S  90°)  «  hSAT  .  (Alt) 

where  the  standard  reference  ellipsoid  is  used  to  approximate  the  sea-level  surface 

required  for  the  calculation  of  the  geocentric  distance  of  the  sub  point. 

The  latitude  of  the  probing  point  is  given  as  ifi..  «  are  sin  <C )  ,IQ).  where  Q.,  is  the 

4  W  4  4 

north-polar  component  of  Q  in  the  standard  inertial  reference  frame* 

To  illustrate  the  peculuritlcs  of  different  effects  of  elevation  offsets  (t  1  and 
cross -elevation  offsets  (y),  one  can  set  the  normalization  term  (1  -f  «  “  -t  y*)' 
appearing  in  Kq.  (A4)  equal  to  unity,  without  introducing  any  significant  error  in 
the  calculation  of  the  offset  correction. 

Using  the  relation  cos  ■  cos  E  cos  (1 ,  cos  X  *  cos  (K  +  <  )  cos(/)  +  y )  from 
figure  A 2,  this  correction  may  he  obtained  as 

A  cos  x  *  cox  X  -  cos  Xj,  ■■ft  cos  x„  I  tan  E|  +  y  cos  E  lain  (i  I  ,  (A9) 

where  the  upper  (lower)  signs  refer  to  E  >  0  (E  v  0);  \  is  the  uncorreeted  value 
of  x  for  i  m  y  m  0.  and  the  expression 
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sinf)  =  P  j  '  |Sj  X  n  j  |  /  cos  E  =  +  ( 1  -  (cos  *  /  cos  E)  ) 


2.1/2 


with  the  same  assignment  of  signs  (upper  for  E  >  0),  defines  the  quantity  p  as  the 
angle  between  the  geocentric  satellite-position  vector  and  the  projection  of  the  solar 
vector  onto  the  plane  of  the  satellite  orbit  (see  Figure  A2), 

The  relations  |Eq.  (A9)]  and  |Eq.  ( A 1 0) J  show  that  for  each  orbit  there  are  four 
characteristic  positions  distinguished  by  vanishing  or  extreme  values  of  the  offset- 
related  corrections  as  follows. 

For  one  pair  of  opposite  positions  defined  by  sin  0  j  =  sin  p  ^  -  0,  only  the 
t  -offsets  have  an  effect  on  \.  For  these  two  positions  the  zenith  angle  goes  through 
its  extremes,  that  is.  xo  j  =  Xq  min  =  |e|  (generally  not  important  since 
hQ  =  hSAT*  and  2  =  xo  max  =  lB0°  ~  lEl  (satellite  either  in  the  dark  or  at  posi¬ 

tion  for  lowest  positive  value  of  h^  near  sunset  or  sunrise  depending  on  specific 
values  of  h.^.j,  and  E).  For  these  positions,  the  correction  equals 
A  cos  \  =  -t  cos  tan  E  (see  Figure  A3). 

The  other  pair  of  characteristic  positions  is  defined  by  the  condition 
sin  p  ,  .  =  ±1.  For  these  positions,  only  the  y -offsets  have  an  effect.  It  is  given 
as  A  cos  X  =  ±  y  (since  XQ  3  =  X,,  4  =  90”)  (see  Figure  A4). 

From  the  foregoing  relations  for  the  angle  p  and  the  offset  coordinate  axes 

-*  a 

t  j  and  y  Jt  one  finds  that  the  expression 


cos  (?)  3  ±  I  tan  E|  tan 


-  90°) 


(All) 


defines  the  quantity  ©  as  the  angle  between  the  "vertical  entrance  plane"  of  the 
instrument  (containing  the  vectors  t.  and  Q)  and  the  "meridional  entrance  plane" 

-*  -4  * 

(containing  the  vectors  tj  and  t  j).  This  meridional  entrance  plane  also  represents 
the  "plane  of  dispersion"  of  the  gratings  in  all  EUV'S  monochromators.  The  entrance 
plane  formed  by  tj  and  y  ^  is  called  "sagittel  entrance  plane"  (parallel  to  the  ruled 
lines  of  the  planar  diffraction  gratings). 

The  Eqs.  (A9)  and  (All),  combined  with  the  relation  A  cos  \  =  -  sin  Ax. 
can  be  used  to  express  the  "vertical  angular  shift"  due  to  offsets  of  e  /  0  and 
p  4  0  in  a  more  comprehensive  form,  that  is. 


Ax  =  -  t  cos  ©  -  y  sin  @  , 

where  cos  ©  was  given  by  Eq.  (All)  and  sin  ©  is,  correspondingly,  given  as 
sin  W  =  ±  1 1  -  (cos  Xq/cos  E)2] 1  ^2  /  sin  xo  * 


(A12) 


(A13) 
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The  corresponding  effect  in  the  re  mote -probing  height  is  then  given  as 

Ahy  •  (e  cos  4  +  y  sin  H)  P  sin  U  ,  (A  14) 

where  U  *  x  ~  90°  is  the  nominal  solar  depression  angle,  the  quantity  d  =  P  sin  U 

represents  the  distance  of  the  remote-probing  point  from  the  satellite  instrument, 

and  W  is  the  angle  between  the  vector  Q  and  the  «  -offset  coordinate  axis.  Since 

this  angle  is  not  required  to  any  accuracy  better  than  a  few  degrees,  the  simplified 

calculation  of  its  value  (for  X  c  X0  and  a ^  «  ±  n^)  implied  by  the  expressions  for 

cosW|Kq.  (All)]  and  sin  W  | Kq.  (A  13))  above,  is  perfectly  adequate.  Our  equation 

for  Ah  .  | Kq.  (A  14)]  is  equivalent  to  Kq.  (6)  of  the  description  of  the  EUVS  experi- 
V 

ment  published  before  launch  [llinteregger  et  al  (1913)  noting  that  the  Kq.  (5)  in  the 
same  report  should  show  cos  0  instead  of  tan  0  written  erroneously). 

This  section  is  concluded  by  giving  an  illustrative  example  in  which  the  required 
corrections  for  offset  pointing  approach  a  magnitude  characteristic  of  "worst  cases,  " 
that  is,  AE-C  observations  made  from  high  satellite  altitudes  and  with  solar  de¬ 
pression  angles  close  to  30°.  Taking  \  =  118°,  P  =  8500  km,  and  K  =  60”  with 

t  23  =  -0.  00103,  y2,j  =  0.  00022  for  MN  No.  23  and  t  (2  =  -0,00002,  >  [2  =  0.  00124 
for  MN  No.  12,  respectively,  one  finds,  with  4  =  22.  935°,  and  d  =  3990  km) 

Ah^  (MN  No.  23)  =  -5.  72  km  and  Ah^  (MN  No.  12)  i  +  1.  85  km  corresponding  to  a 
difference  of  almost  8  km  for  the  corrected  probing  heights  for  simultaneous 
observations  by  MN  No.  23  and  MN  No,  12. 

In  Figure  Al,  the  plane  of  drawing  (a)  is  perpendicular  to  the  direction  toward 
the  center  of  the  solar  disk.  The  offset-pointing  subsystem  allows  selection  of  any 
one  of  256  target  points.  The  elevation  and  cross  elevation  offset  angles  relative 
to  central  solar  pointing  are  given  as  t =  (NK-8)At  and  Fgpg  =  (NG-H)Ay, 
respectively  (At  =  Ay  =  0.001018  rad/step).  The  positive  directions  of  the  <  and 
y  axes  are  uniquely  determined  from  the  known  directions  of  the  spacecraft  spin 
axis  (unit  vector  of  +  ZC„  axis  equal  to  outward  pointing  unit  vector  of  SPS  azimuth 

4  4 

axis  designated  a.)  and  unit  vector  of  direction  to  center  of  solar  disk  (s.),  for 
example,  Vj  =  Sj  X  aj/sinG,  as  shown  in  the  schematic  (b).  Since  the  spin  axis 
a  t  was  routinely  held  within  ±  1-2°  of  the  normal  direction  to  the  orbit  plane,  one 
finds  G  ^  90°  -  |e|  with  E  >  0  (E  <01  expressing  the  solar  elevation  over  the 
satellite's  orbit  plane  (E  <0  when  sun  is  on  the  negative  side  of  that  plane),  that  is, 

*4  4  4  .  *4 

setting  aj  =  n(  sign  (^  •  s().  where  n ^  is  the  normal  vector  of  the  orbit  plane. 
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Various  quantities  used  in  the  text  are  illustrated  in  Figure  A  2.  SC  =  space¬ 
craft.  O  -  center  of  earth.  Q  =  remote  probing  point.  R  =  radius  of  earth, 
hQ  =  height  of  remote  probing  point.  \  =  solar  zenith  angle  of  the  target  point. 

\0  solar  zenith  angle  of  the  center  of  the  solar  disk. 

Figure  A3  shows  the  spacecraft  (SC)  at  the  position  where  0  =  0°.  \  =  180-E 
Thus  the  y -offsets  have  no  effect  on  \. 

Figure  A4  shows  the  spacecraft  (SC)  at  the  position  where  0  =  90°.  \  =  90°. 

Thus  the  t  -offsets  have  no  effect  on  \  . 
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Figure  A2.  Various  Quantities  Used  in  Text 
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A 2.  EFFECTS  OF  FIKLDS-OF-V IEW 


The  effects  of  the  field-of-view  (FOV)  for  the  so-called  "full-disk  mono¬ 
chromators"  and  "reduced-field  monochromators"  are  discussed  separately  below. 
The  monochromators  of  the  first  group  had  been  designed  with  emphasis  on  the 
measurement  of  incident  EUV  fluxes  at  all  wavelengths,  that  is,  using  a  wavelength¬ 
scanning  mode  of  128-steps  for  each  of  the  scan-capable  monochromators.  How¬ 
ever,  those  scan  observations  have  been  scheduled  only  for  those  parts  of  an  orbit 
for  which  the  combination  of  solar  zenith  angle  and  satellite  height  has  been  favor¬ 
able  in  the  sense  of  expecting  only  reasonably  small  corrections  for  atmospheric 
attenuation.  For  the  remaining  time  of  instrument  operation,  all  monochromators 
have  been  kept  in  fixed  wavelength  settings,  corresponding  to  certain  standard  steps 
(step  64  for  all  EUVS-1  MN's.  steps  64,  24.  88,  S>6.  or  124  for  all  EUVS-2  MN's) 
to  provide  the  data  desired  for  EUV  absorption  analysis  of  atmospheric  structure. 
Conversely,  the  reduced-field  monochromators  had  been  designed  with  primary 
emphasis  on  EUV  absorption  analysis,  and  their  use  for  observing  temporal  varia¬ 
tions  of  incident  solar  EUV  at  these  fixed  wavelengths  reflects  only  a  secondary 
objective  in  terms  of  aeronomical  interest.  The  latter  interest  would  of  course 
have  preferred  full-disk  flux  data  rather  than  observations  restricted  to  some  small 
central  part  of  the  area  of  EUV  sources  on  the  solar  disk. 

A2. 1  Full-Dink  Monochromator* 

Approximately  one-half  of  the  monochromators  used  in  the  EUVS  instrument  on 
the  Atmosphere  Explorer  satellites  AE-C,  D,  and  E  are  called  "full-disk  MN's." 
since  their  unencumbered  field-of-view  extends  beyond  the  solar  disk  in  both  the 
t  -direction  (unencumbered  meridional  acceptance  width  of  Ar  ^  ‘  0.  0174  rad)  and 
the  y  -direction  (sagittal  width  of  Ay^  i  0.0174  rad).  As  can  be  seen  in  Figure  Al, 
those  unencumbered  widths  are  greater  than  the  angle  subtended  by  the  photospheric 
disk  (Ae,,^  =  AyslJ^  i  0.0002  rad).  This  excess  of  instrumental  acceptance 
widths  is  sufficient  to  guarantee  full-disk  response  even  for  observations  made  with 
intended  or  inadvertent  offset  angles  greater  than  those  actually  used  in  all  geo¬ 
physical  routine  turnons  of  the  EUVS  experiment  on  all  of  these  AE  satellites.  The 
center  of  the  instrumentally  viewed  field  of  solar  EUV  source  is  then  of  course 
identical  with  the  center  of  the  solar  disk,  and  the  algorithms  developed  in  Section 
Al  are  relevant  only  for  the  assessment  of  the  spread  of  zenith  angles  (A\)  and 
remotely  probing  heights  (Ah^)  associated  with  the  spread  of  solar  source-point 
coordinates,  0  v  |c  |  v  0.0046  and  0  *■  |y|  ^  0.0046,  regardless  of  SUS  target 
offsets  if  any  (up  to  at  least  2-8  offset  steps  as  illustrated  in  Figure  Al)  and  regard¬ 
less  of  individual  monochromator  misalignment  angles  of  values  even  exceeding  the 
worst  actual  cases  of  t  4  0  and  y  0  in  the  AE-C  instrument. 


For  those  full-disk  MN's  the  expression  (Eq.  Al4)  for 
replaced  by  the  simple  relation 


is,  therefore, 


(A15) 


where  P  cos  (180°  -  \  )  is  the  distance  of  the  remote-probing  point  from  the  satel 
lite,  and  the  total  spread  of  zenith  angles  associated  with  all  source  points  on  the 
solar  disk  becomes 


(A16) 


for  all  values  of  the  solar  zenith  angle  associated  with  the  ray  from  the  center  of 
the  solar  disk  )• 


